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SUMMARY 


The influence of changes in the precursor region flow properties (resulting 
from absorption of the radiation from the shock layer) on the entire shock-layer 
flow phenomena around a Jovian entzry body is investigated under physically 
realistic conditions. In the precursor region, the flow is considered to be 
inviscid and the variations in flow properties are determined by enploying the 
small perturbation technique as well as the thin- layer approximation. The 
flow in the shock layer is assiamed to be steady, axisymmetric , and viscous. 

The analysis is carried out by considering both the chemical equilibrium and 
nonequilibrium composition of the shock -layer gas. The effects of transitional 
range behavior (slip boundary conditions on the body surface and at the shock 
wave) are included in the analysis of high altitude entry conditions. 

Realistic thermo-physical and radiation models are used, and results are 
obtained by employing the implicit finite difference technique in the shock 
layer and an iterative procedure for the entire shock layer-precursor zone. 
Results obtained for a 45° hyperboloid blunt body entering Jupiter's atmosphere 
at zero angle of attack indicate that preheating of the gas significantly 
increases the static pressure and temperature ahead of the shock for entry 
velocities exceeding 36 Jon/sec. The nonequilibrium radiative heating rate 
to the body is found to be significantly higher than the corresponding equi- 
libritam heating. The precursor heating, in general, increases the radiative 
and convective heating to the body, and this increase is slightly higher for 
the nonequilibrium conditions. 
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1. INTRODUCTION 


A space vehicle entering a planetaury atmosphere encounters a wide 
range of flow conditions ranging from free molecular flow at high 
altitudes to continuum flow at low altitudes. Since experimental 
facilities caumot adequately simulate conditions expected during entry 
into the outer planetary atmospheres, most of the required information 
must be obtained from theoretical studies. 

During the high speed entry, the atmospheric friction works as a 
brake to slow the spacecraft and the gas around the body in the formed 
shock layer is heated by the dissipated kinetic energy. Radiation plays 
a very importemt role in the euialysis of flow phenomena around an entry 
body at high speeds. In many instances, the radiative energy trans- 
ferred to the body from the high temperature shock layer gas exceeds the 
convective euid aerodynamic heat transfer. Radiative energy transfer 
from the shock layer of a blunt body into the free stream reduces the 
total enthalpy of the shock layer while increasing the enthalpy of the 
free stream gases. Because of this increase in enthalpy, the entire flow 
field ahead of the shock layer and around the body is influenced signi- 
ficantly. The phenomena of chsuige in flow properties ahead of the 
shock wave due to the energy interaction from the shock layer is called 
the "praecursor" or "praecurrere" (prae = before + currere = run) which 
mecuis "forertinner. ” In the present context, therefore, the precursor 
flow region is considered to be the region cdiead of a shock wave in 
which the flow field parameters have been changed from free stream condi- 
tions due to absorption of radiation from the incandescent shock layer. 
Most of the radiative energy transferred from the shock layer into the 
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cold region ahead of the shock is lost to infinity unless it is equal to 
or greater than the energy required for dissociation of the cold gas. 
tJhen the photon energy is greater than the dissociation energy, it is 
strongly absorbed by the cold gas in the ultraviolet continuum range. 

The absorbed energy dissociates and ionizes the gas and this results in 
a change of flow properties in the precursor region. In particular, the 
temperature and pressure of the gas is increased while velocity is 
decreased. The chcuige in flow properties of the precursor region, in 
turn, influences the flow characteristics within the shock layer itself. 
The problem, therefore, becomes a coupled one and iterative methods are 
required for its solution. 

Only a limited number of analyses on radiation induced precursor 
flow is available in the literature. Works available vintil 1968 are 
discussed, in detail, by Smith [1,2]*. By employing the linearized 
theory of aerodynamics. Smith jLnvestigated the flow in the precursor 
region of a reentry body in the earth ' s atmosphere . The cases of plane , 
spherical, and cylindrical point sources were considered and solutions 
were obtained for a range of altitudes and free stream conditions. It 
was found that for velocities exceeding 18 km/sec, precursor flow effects 
are greatest at altitudes between 30 and 46 km. It was further concluded 
that preheating of air may cause an order of magnitude increase in the 
static pressure and temperature ahead of the shock wave for velocities 
exceeding 15 lon/sec. A few other works, related to the effects of up- 
stream absorption by air on the shock layer radiation, are discussed by 
Liu [3,4] . Some works on precursor ionization for air as well as 
hydrogen-helium atmosphere are presented in [5-9] . 

* The numbers in brackets indicate references. 
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In the analysis of most shock— layer flow phenomena, the contribution 
of radiation-induced precursor effects usually is neglected. Garrett [10] 
presented a detailed review of the various methods used for solving the 
radiating flow field at the stagnation region. Also, various methods of 
solution of radiating shock layer are discussed by Anderson [11] . 

Sutton [12] separated the radiating flow field into an outer inviscid 
layer and an inner boimdaary layer; the two solutions are coupled by 
radiative transport through both layers and by the ' boundary displace- 
ment thickness. Kumar, Tiwari, and Graves [13] considered the entire 
shock layer as viscous flow region and used a time -dependent method to 
obtain the solutions for small angle of attack. Davis [14] presented a 
method for solving the viscous shock-layer equations for stagnation amd 
down stream flow. Moss [15-17] applied this method of solution to 
reacting multicomponent mixtures. The precursor effects were neglected 
in all the above studies. However, a limited number of studies which 
include this effect are available in the literature. Lasher and 
Wilson [18,19] investigated the level of precursor absorption and its 
resultant effect on surface radiation heating for e£u:th*s entry condi- 
tions. They concluded that, for velocities less than 18 km/sec, pre- 
cursor heating effects are relatively unimportant in determining the 
radiative flux reaching the sxarface. At velocities greater than 
18 km/sec, the amount of energy loss from the shock layer and resultant 
precursor heating correction was found to be significantly large. 

Liu [3,4] also investigated the influence of upstream absorption by 
cold air' on the stagnation region shock layer radiation. The thin layer 
approximation was applied to both the shock layer and the preheating 
zone (the precursor region) , The problem was formulated for inviscid 
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flow over smooth blunt bodies but the detailed calculations were carried 


out only for the stagnation region. The general results were compared 
with results of two approximate formulations. The first approximate 
formulation neglects the upstream influence and the second one essen- 
tially uses the iterative procedure described by Lasher and Wilson 
[18,19]. The results eire compared for different values of a radiation/ 
convection parameter. 

As mentioned eeirlier, the cold gas absorbs energy only by photodis- 
sociation and photoionization in the precursor region. The absoi:ption 
coefficients are a continuous nonzero function of photon energy (because 
of bound-free transition) for all values of photon energy exceeding the 
dissociation potential of the molecule. A critical review of ultra- 
violet photoabsorption cross sections for molecules of astrophysical cind 
aerohomical interest, available in the literature up to 1971, are given 
by Hudson [20] . Specific information on photoionization and absorption 
coefficients of molecular hydrogen is available in [20,21] . 

In the shock layer region, the gas may be treated as gray or non- 
gray, Anderson [11] concluded that a gray gas analysis is not suffic- 
iently accurate for entry applications and suggested use of nongray 
models. The frequency dependence of the absorption of coefficient for 
a nongray gas may be treated either in detail or by a "step model." 

There exists several ccmputer programs for the detailed fretjuency 
dependence of the absorption coefficient which are developed by Nicolet 
[22], Wilson [23], and Thomas [24], In a step model, the frequency 
dependence is broken into a number of discrete steps, Falanga and 
Olstad [25] presented a 38-step model for 90% CO^ and 10% (by 
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voliame) mixture which included 15 steps to model the continuum and 23 
steps to model the line contribution to the radiation transport, Zoby, 
Sutton, and Moss [26] developed a 58-step model for hydrogen and heliimi 
mixture. The transitions considered in this model are: the bound-bound, 

bound-free and free-free transitions for atomic hydrogen, the boxind-free 
and free- free transitions for the negative hydrogen ion, and the Lymsin ■ 
and Werner band systems for molecular hydrogen. This 58-step model is 
fairly accurate and compares very well with the results of Nicolet's 
detailed model for hydrogen/helium species [27] , 

The total radiative transport is an integral over both the fre- 
quency spectrum and the physical space. The methods for calculating the 
divergence of the radiative flux and other conservation equations are 
available in [28-30] , 

It is very well documented in the literature that the degree of 
rarefraction of a flow is measured by the Reynolds number. Therefore, 
for a fixed blunt body at low altitudes where the Reynolds number is 
high, the flow will behave like an ordinary viscous flow which lies 
within the scope of the Navier-Stokes equations. At higher altitudes, 
where the Reynolds number is low, the theory of free molecular flow can 
be used. The transition zone between these two regions has been divided 
into several stibregions which are discussed in greater detail by Hayes 
and Probstein [31] , and Cheng [32] has provided the different methods of 
solution valid within each region. In the continuum range, the flow 
phenomena is investigated through use of the Navier-Stokes equations. 

In the transition reinge (from the continuum end), however, use of the 
Navier-Stokes equations is still justified for the main flow field but 
the boxmdary conditions cannot be satisfied in the usual manner. Thus, 
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the characteristic feature of flow of a slightly rarefied gas, which 
sharply distinguishes it from the continuxim flow, is the chcuage in the 
boundary conditions at the body surface [33] and shock wave | [^4-35] . 
Instead of using the Rankine-Hugoniot conditions as boundary conditions 
at the shock wave, Probstein 2 uid Pam [34,35] introduced the concept of 
••shock wave slip" aai interpretation of the trzmsported effects behind 
the shock. Rott emd Lenard [36] have shown that the effects of velocity 
slip and temperature jump on the body surface cannot be neglected in 
compeurison with other low Reynolds number corrections. A semi-macro- 
scopic argument which leads to the simple expression for velocity slip 
and temperature jump is given in [33,36] . 

From the literature survjey, it is quite clear that no work is avail- 
ed>le which considers the influence of precursor heating on the shock 
layer flow phenomena around a Jovian entry body. A few studies that are 
available deal only with the case of chemical equilibrium in the shock 
layer; the case of shock layer chemical nonequilibrium flow has not been 
considered. Also, no' consideration has been given to investigate the 
transitional range shodk layer flow phenomena encountered at high Jovian 
entry altitudes. In an actual entry situation, the influence of pre- 
cursor heating, nonequilibriimi chemistry in the shock layer, and transi- 
tional range flow phenomena may be strongly coupled. Thus, it is 
essential to investigate the extent of influence of each phencxnena 
separately and jointly in order to assess the true behavior of flow 
aucound the entry body. This information is of vital importance in 
determining the convective and radiative heating of the entry body. 

The main purpose of this study, therefore, is to investigate the 
influence of changes in the precursor region flow properties on the 
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entire shock layer flow phenomena aroimd a Jovian entry body. The cases 
of shock layer chemical equilibrium as well as chemical nonequilibrium 
are considered, and the effects of transitional range behavior are 
included in the analyses of high altitude entry conditions. In order 
to accomplish these objectives in a systematic manner, the entire problem 
has been divided into four subproblems as: 

1. Investigation of the radiation induced precursor region flow 
phenomena . 

2. Effects of shock and body slip conditions on viscous equilibriim 
flow. 

3. Influence of precursor heating on viscous equilibriijm flow. 

4. Influence of precursor heating on viscous nonequilibrium flow. 

Basic foimiulation of the entire problem is presented in Chap. 2. 

The boundary conditions for different flow regimes are given in Chap. 3. 
Information on thermodynamic and trem sport properties for each species 
considered in different flow regimes are given in Chap. 4, Information 
on chemical reactions and reaction rates for both ec[uilibrium and non- 
equilibrium conditions are given in Chap, 5. Discussions on radiation 
models and radiative flux equations are presented in Chap. 6. Solution 
procedures for the precursor and shock layer regions are discussed, in 
some detail, in Chap. 7. Discussions of all results are presented in 
Chap , 8 . 
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2 . BASIC FOBHULATION 


The physical model and coordinate system for a Jovieui entry body is 
shown in Fig. 1. The entire flow field adiead of the body c«ui be divided 
essentially into three regions: the free stream, the precursor region, 

and the shock layer. The flow properties are considered to be uniform 
at large distauices from the body. In this section, governing equations 
are presented for the precursor as well as shock layer region. However, 
it would be appropriate here to discuss first the Joviam atmospheric and 
entry conditions. 

2.1 Free-Stream Region 

Information on Jupiter's atmospheric conditions are availaUble in 
[37-39] . In the past, the ncminal composition of the atmosphere was 
assumed to be 85 percent hydrogen and 15 percent helium by mole fraction. 
Recently, this has been changed to 89 percent hydrogen and 11 percent 
heliimt [39] , For different altitudes of entry, the free-stream condi- 
tions used in this study are given in Tables 1 and 2. The temperature 
of the atmosphere (i.e., T^) is taken to be constant at 145 K and the 
free-stream enthalpy ceui be calculated by following the procedure given 
in [1,2] as 

H = 1.527 R T <2.1) 

GO 00 

where R = 8.315 Joules/ K-mole is the universal gas constant. The 
number density of hydrogen can be calculated by the ideal gas law and 
the relation caui be expressed as 

N = (7.2431172 X 10^^) (P /T )X„ (2.2) 

“2 2 

where is the mole fraction of and P^ has the units of N/ia^. 
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Table 1 Altitude and free-streain Jupiter entry conditions 
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Table 2 Free-sbream and shock conditions for Jovian entry. 


Free stream 

V^, km/sec 

*s' 

q(0), erg/cm^ 

Z = 95, km 

38 

16,610 

1.35 E12 

p - 1.29 E-3, kg/m^ 

09 

35 

15,400 

7.75 Ell 

P = 673, N/m 

32 

14,080 

3.52 Ell 


30 

13,550 

2.01 Ell 

Z * 103 

40 

16,890 

1.16 Ell 

p =8.56 E-4 
00 

35 

15,040 

4.70 Ell 

P = 448 
00 

33 

14,250 

3.28 Ell 


30 

12,810 

1.142 Ell 

Z = 116 

45 

18,227 

1.09 E12 

p =4.65 E-4 
^00 

39.09 

15,886 

4.76 Ell 

P » 244 
00 

35 

14,480 

2.18 Ell 


30 

12,480 

4.87 ElO 

Z = 131 

43.21 

16,390 

3.86 Ell . 

p = 2.32 E-4 
00 

38 

15,210 

1.61 Ell 

P = 122 
00 

35 

13,880 

8.72 ElO 


30 

12,030 

1.90 ElO 

Z = 150 

42 

15 , 050 

9.60 ElO 

p = 9.29 E-5 
^00 

40 

14,520 

6.96 ElO 

P = 49 

GO 

35 

13,140 

2.57 ElO 


30 

11,600 

6.20 E9 
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2.2 Precursor Region 

In this region, the flow is considered to be steady eind inviscid. 
To investigate the changes in flow properties in this region, both the 
small perturbation theory of classical aerodynamics eind the thin~layer 
approximation of hypersonic flow have been used in this study. Funda- 
mental principles of these approximations are briefly discussed in the 
following subsections. 

2.2.1 Small Perturbation Theory 

For application of the small perturbation theory, basic conserva- 
tion equations for the precursor region can be written as [40,41] 


Hass Continuity: 

V • CpV) =0 (2.3) 

Homenttim : 

p(V • VV) = -Vp (2.4) 

Energy: 

p(V • VH^) = (2.5) 

Species Continuity: 

p (V • VC ) » K (2.6) 

a a 

State : 

p = pRTZ CC /W ) (2.7) 

^ a a a 


where the total enthalpy per unit mass is given by 

H » H + v2/2 
T 

In the above equations, Q = '^•q is the net rate of radiant energy 

R R 

cibsorbed per unit volume per imit time, represents the net rate of 
production of species a per unit volume per unit time, and w^ is the 
raole<^lar weight of species a. 
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AS a result of Increased fluid enthalpy, the entire flow field in 
the precursor region is perturbed. By following the small perturbation 
technique of classical aerodynamics, the flow properties can be expressed 
in perturbation series as [1^2, 40-43] 


p = 

p, tt + + 2 + ...) 

(2.8a) 

p = 

P. (1 + + Pj + ...) 

(2.8b) 

V = 

+ ''l + '^2 ' 

(2.8c) 

H = 

H + CH, + H, + ...» 

OO OO ^ 

(2.8d) 

T = 

+ Ti + T 2 + ... . 

(2.8e) 

C = 
a 

C + C ^ + C ^ + ... 

a« al a2 

(2.8f) 


In these equations, all the perturbation vcuriables Cexcept temperature) 
are expressed in nondimensional form. The unit vector k represents the 
direction of unperturbed free-stream velocity. 

If and can be considered as first-order perturbation terms, 
then substitution of Eqs. (2.8) into Eqs. (2.3)- (2.7) results in the 
first-order perturbation equations as 
Continuity • 

7 • Vi + 3p^/3z = 0 (2.9) 

Monentum : 

3V^/3z = -(1 /ym2)Vp (2.10) 

Energy : 

3H /3z = Qj^/(p„ v3) (2. IX) 

Species: 

3C /3z = K / (p V ) (2.12) 

a 00 00 

where 


= H- + V, 

1 Iz 


and Y represents the ratio of specific heats. 


(2.13) 
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The bOLindary conditions eure that perturbation quantities veuiish at 
z 00 and that no singularities exist except at the origin. 

The radiation effect on the gas ahead of the shock produces H, 
and electrons e by photodissociation emd photoionization, and also 
increases the enthalpy. Any other species which may be produced aire 
neglected. The contribution of radiation to the gas pressure is neglected. 
It is further assumed that the internal degrees of freedom of various 
species (i.e., vibrational and electifonic modes) are not exited. For 
this gas model, the ec[uation of state (for the first order perturbation) 
can be expressed as [1,2] 

P. * (400/180.17) [(C + C )/2] + (T./T) + p. (2.14) 

J. n "2'^ ^ ^ 

By following the procedure described by Smith [1,2] , the first-order 
perturbation relation for enthalpy is found to be 

H, « (l/V^) {1.527 RT, + [(5/4)RT + 1/2] C 

1 * 1 * H 2 + 


+ [C3/4)RT + D]C } (2.15) 

H 

where I «md D represent the ionization and dissociation energy res- 
pectively. It should be pointed out here that D in the above equation 
actually represents half the energy required for dissociation. 

As pointed out earlier, the upstream gas absorbs the energy radiated 
from the shock layer in the ultraviolet continuum range. The radiation 
from the perturbed gas due to recombination Ci.e., emission) is neglected. 
The amount of radiative energy cdasorbed by the perturbed gas per unit 
volume and time, Q , is given by 

Q = N„ /“ H a(v)dv (2.16) 

*^R H_ o V 
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where is the number density of is specific irradiance amd 

"2 

a(v) is the photon absorption cross section of at frequency v. 

In determining the rate of production of species in the precursor 
region, only photodissociation and photoionization are considered.. 
Reccmbination is assumed to be a second order effect and, therefore, is 
neglected in the present linearized treatment. The net rate of produc- 
tion of species, therefore, is given by [1, 28] 

-CT /hv) a Cv)dv (2.17a) 

H 1 0 V O 

^ j (H /hv) a (v)dv (2.17b) 

where m represents the weight of an H molecule (in grams per molecule), 

and a (v) and a (v) are the absorption cross section for photodissocia- 
D I 

tion and photoionization, respectively. 

2.2.2 Thin Layer Approximation 

The concept of thin shock layer theory (usually applied to hyper- 
sonic shock layer flows [31] ) is also applied to investigate the precur- 
sor effects. The curvilinear orthogonal coordinate system, shown in 
Fig, 1 is selected and the differential equations for a hypersonic plane 
or axisymmetric flow can be written in the present coordinate system as 


[42] 

C8/3s)Cpur^) + 0/3n)(pvXr^) = 0 (2.18) 

p[uC3u/9s) + (Xv(9u/3n) - Kuv] + (3p/9s) = 0 (2.19) 

p[u(9v/9s) + Xv(9v/9n) + Ku^] + XC3p/3n) = 0 (2.20) 

p [ (u/X) (9H/3s) + vC3H/3n)] + (Xr^ ) “^ [ (3/3n) (Xr^q^^) ] = 0 (2.21) 

p[(uA)OC^9s) + v(9C^/3n) * = 0 (2.22) 


where K = K(s) = 1/R , X =■ 1 + Kn, and j = 0, for plane flows and 1 for 

s 

axisymmetric flows. It should be noted that, according to the notations 
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used in Fig. 1, all quantities appearing in the edxsve equations should 
have a prime superscript (i.e., u',v', p', H*, etc.), 2 uid all physical 
coordinates should have a superscript * (i.e., s*^ n*, r*, etc.). How- 
ever, for the sake of clarity, these notations have been emitted from the 
equations. 

If the precursor region is assumed thin, then one can meJee the 

approximations that (n/R )«1, 3/3s« 3/3n, and r^ is not a function of 

s 

n. In this case X « 1, and Eqs. (2.18)- (2.22) reduce to simpler forms 


as [42] 

0/3n) (pv) - 0 (2.23) 

pvOu/9n) * 0 (2.24) 

pv(3v/3n) + (3p/3n) » 0 (2.25) 

pv(3H/3n) + (3qj^/3n) - 0 (2.26) 

pv(3C /3n) - K « 0 (2.27) 

a a 


The similarity between these equations etnd the small perturbation 
Eqs. (2. 9)- (2.12) should be noted. 

In present application to the hydrogen-helium atmosphere, Eq. (2.27) 
Will be written for atomic hydrogen 2 uid hydrogen ions. In Eq. (2,26), H 
represents the total enthalpy amd is given by the relation 

H - - h + (u2 + v2)/2 (2.28) 

where 

h" * 1.527 RT + [(5/4) RT + I/2]C„+ + [(3/4)RT + D)C„ (2.29) 

“2 

Note that Eq. (2.29) is slightly different thaui the relation for per- 
turbation enthalpy given by Eq. (2.15). 
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2.3 Shock- Layer Region 


In this region, the flow conditions for which the present analysis 
is carried out are axisyinmetric , steady, laminar, viscous, and compress- 
ible. It is fiirther assumed that the gas is in local thermodynamic 
equilibrium auid the tangent slab approximation is valid for radiative 
transport. The reacting multicomponent gas mixture is treated in both 
chemical equilibrixm and nonequilibrium conditions. 

2.3.1 Chemical Equilibriim 

The viscous shock layer conservation equations presented in [14-17] 
are 'a set of equations that are valid uniformly throughout the shock- 
layer region. The methods of obtaining these equations are discussed 
in detail in those references. First the conservation equations are 
written for both the inviscid cind the boundary- layer regions in the body- 
oriented coordinate system. Then these equations are nondimensionalized 
in each of the two flow regions with variables which are of order one. 
Terms in the resulting sets of equations are retained up to second order 
in the inverse square root of Reynolds number. Upon combining these two 
sets of equations, so that terms up to second order in both regions are 
retained, a set of equations uniformly valid to second order in the 
entire shock layer is obtained. The nondimensional form of the viscous 
shock-layer equations that are applicable in the present case can be 
written as 

Continuity: 

(3/3x) [ (r+y cos 0)pu] + (3/3y) [(l+yK)(r+y x cos,9)pv]=0 (2,30) 

X-momentum r 


/ u 3u 3u uv< \ 1 

\1 + yic 3x ^ ^ 3y ^ 1 + yK/ 1 + j 

“ ^ - rT7^)J " n rr 

./ii _ } 

\3y 1 + yK / 


cos 6 

r+y cos 0 


(2.31) 
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Y-mGmen‘tum: 
P 


r — ^ ^ + V ^ ^ - 0 

\1 + y< 3x 3y 1 + y / 3y 


(2.32) 


Energy : 


where 


u ^ 

y)c 3x 


3H ^ 3H\ 

3x 3y ) “ 


^ . ecu; 


3y 1 + yic 
3C. N 


If t 


I 3y jpr 3y * Pr i 3y " ‘’i P 

N 

" I / K COS 9 \ r u 3H _u_ j,s , 

1 + yic J \ 1 + yV r + y cos 0 / jpr ?y ” Pr “i 3y 

h. J. +-H. (pr-l)u |a - J£^] } .fi!« 

^ ^ Pr 3y 1 + yKj * 3y 

+ q + H 2 £ J .')1 

^R\l + yjc r J \ 


(2.33) 


H = h + uV2. 

The terms used to nondimens ionalize the ahove equations are defined as 


X =» x*/R* 
n 

V » v*/V* 

' 09 

Pr = C*p*/K* 
P 

y * y*/R* 
n 

p » pVp* 
00 

Le. . * p*C*D*./K* 
13 P 13 

r = r*/R* 
n 

p = 

L. . = p*C*or ,/K* 
13 P 13 

K = KV(U^^jC*„) 

C = C*/C* 
p p P" 

P = P*/(P^*^) 


u . = P*(V*2/C* ) 
ref w ' poo 

T = T*C /V*2 

poo' GO 

h = h*/V*2 
00 

K *= K*/R* 
n 

u * u*/v* 
00 

J. - j*R*/y* ^ 
1 in ref 


|l/2 
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In Eq. (2.33), represents the mass fltix relative to the mass 


average velocity amd is given by the expression [14,44] 

[ NI ■ -I 

I (3Cj^3Y) + (L^/T) (3T/3Y)J 


(2.35a) 


^ere 


iK 


Le. , i » 


NI C . NI 


Le « Z ( -i) Z CC./M.L. .) 
i .AM./._ t 3 XT 

j*l 3 3=1 J J J 


D?^l 


j?^l 


Ab 


iK 


i«. -{ 


(M./M)Le.„ + [1-(M./M„)1 

X Xiv X Jx 


NI 

z 

j=i 


Le. .C. I 
13 3 / 


The last term in Eq. (2.35a) represents the contribution of thermal 
diffusion. The quamtity Le^^ represents the multi-component Lewis 
number, and represents the binary Lewis Semenov numbers; both are 

defined in Eq. (2.34). If thermal diffision cam be neglected and L^^ can 
be taken as constant for all species, then Eq. (2.35a) reduces to 

J. = - (v/Pr)L. . C3C./3y) (2.35b) 

In the present study, use is made of Eq. (2.35b), and the value for 
L^^ is taken to be 1.1 [45,46]. 

The ejq>ression for the equation of state for a hydrogen/helium 

mixture is given by Zoby et al [47] as 

T* = C^[ (p*/l0l32 50 )V(P VO. 001292)^] (2.36a) 

H* = C[(pVl013250)“/ (p/0. 001292)”] (RT /M) (2.36b) 

n O 

where 

K = 0.65206 = 0.04407 £n(X .) 

H2 
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I - 0,67389 - 0,04637 An(X„ ) 

“2 

m - 0.95252 - 0.1447 £n(X_ ) 

**2 

n - 0.97556 - 0.16149 ln{X„ ) 

“2 

a - V sin e El + 0.7476 (1-X„ )] 
t - 

CTO « - 545.37 +61.608 - 22459 uj + 0,039922 

- 0.00035148 + 0.0000012361 u| 

CHU « 5.6611 - 0.52661 o£ + 0.020376 - 0.00037861 u£ 

b c c 

+ 0,0000034265 - 0.000000012206 

C_ - CTO + 61.2 (1-X„ ) 

T 

C„ - CHU - 0.3167 a-X„ ) 

H 

and represents the mole fraction of H . 

X2 2 

The set of governing equations presented above has a hyperbolic/ 
parabolic nature. The hyperbolic nature enters through the normal 
momentum equation. If the shock layer is assumed to be thin, then the 
normal momentum equation can be expressed as 

pu^K/a + y<) - Op/3y) (2.37) 

If Eq. (2.32) is replaced with Eq. (2.37), then the resulting set of 
equations is parabolic. These ec[uation8 can, therefore, be solved by 
using numerical procedures similar to those used in solving boundary* 
layer problems [14,15], 

2.3.2 Chemical Nonequilibrium 

For the condition of chemical nonequilibrium, the basic governing 
equations (continuity, X*momentum, Y-momentum, emd energy) are essen* 
tially the same as given for the chemical equilibrium condition. The 
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«p«cl«s continuity oration, howaver, is needed and this is given by 
the relation 

/ u aCiS . 

(rr^ 9^ r \ - 

( 3 [ Cl + nic) (r + n cos 0)^J^1 ^ (2,30) 

(i + nic) (r + n cos 6)^ ^ 

where represents the rate of production of chemical species in the 

shock layer. The equation of state given by Eq. (2.36) is valid only 

for the chemical equilibrium case. For the case of chemical nonequi* 

librium, the equation of state is given by the relation [28) ^ 

P*V* - E (N. ) R*T* (2.39) 

i ^ 

where is the number of moles for the i-th species. This result is 
reminiscent of the thermal equation of state for a perfect gas. The stmi 
in parentheses, however, is not a constant since the total nuadser of 
moles change as the chemical balamce changes. 
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3. BOUNDARY CONDITIONS 

As pointed out earlier, the slip boundary conditions eure not 
importeuit at low altitudes but they cannot be neglected at higher 
altitudes. Since both the slip amd no-slip conditions have been used 
in -this study, they will be discussed sepeirately in this chapter. 

3.1 No-Slip Bound 2 u:y Conditions 

At tile body surface (wall) , no velocity slip zmd no temperature 

jump are assumed. Consequently, the velocities at the surface are 

V - 0 (3.1) 

u - 0 (3.2) 

The wall temperature for this study is specified as 

T = constant (3.3) 

w 

The Rankine-Hugoniot relations cure used to determine the flow 
properties immediately behind the shock. The nondimen sional form of 
the shock relations can be written as [45] 

Continuity i 

p V _ = -sina (3,4) 

s- s 

Momentum: 


u' _ = sina (3,5) 

s 

Ps_ * Pg+ + Sinead- l/pg-) (3.6) 

Energy : 

h _ =* h + + (sin2a/2) (1-1/p ^_) (3.7) 

s s s 

where a is shown in Fig. 1, and u^ and v^ are velocity components 

expressed in a shock-oriented coordinate system. The relations for 

u and V in the body-oriented coordinate system ceui be written as 
s s 


u = u* sin (a + B) + v* cos (a + 3) 
s s s 


(3,8) 
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(3.9) 


V =s - u' cos Cci + 6) + V* sin (a + B) 
s s s 

where angle B is indicated in Fig. 1. 

3.2 Slip Boundeury Conditions 

Shidlovsky [33] has shown that at the body surface the velocity 

slip and temperature jump conditions are of the same order as the Knudsen 

number. The Knudsen number, K^, is defined as the ratio of the particle's 

mean free path I eind the characteristic dimension L of the body (i.e., 

K = A/L) . The ordineury boundary conditions (which correspond to 
n 

continuum conditions) are obtained when -► 0. However, for the 
transitional range (i.e., for -*■ 0(1)), in order to be consistent 
with the Navier-Stokes equations of motion, a linear relation between 
the conditions at the wall and the flow should be assumed. This ccui be 
done by a semi-macroscopic argument which leads to the simple expression 
for velocity slip and temperature jump as [48-51] 


u = 

A^ (y/P) (P/p)^'^^ Ou/3y) 

(3,10) 

T = T 
w 

+ A 2 (K/P (P/p)^^^ OT/3y) 

(3,11) 

V = 0 


(3,12) 


where and are constants and are given by 

Ai - (ir/2)^^^, Aj = (15/8) [(2-02 j/Oj] 

The terms and are slip and thermal accommodation coefficients 
respectively amd are dependent on the nature of the surface and fluid. 
However, in actual flight conditions both and are expected to be 1, 
Since the transport and thickness effects are important at higher 
altitudes, the conditions imposed at the shock cannot be calculated by 
using the classical Rankine-Hugoniot relations, Probstein euid Pan [34] 
have shown that the thickness effect is of a higher order in a and. 
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therefore, .it is neglected in the present study. The information on 

thickness effect can be found in [52,53]. In the present case, the 

shock may still be considered "thin" when compared to the thickness of 

the viscous shock layer. As such, the thin-layer approximation 

C3/3X* )<< 0/3y ' ) and y*/R can be used in the shock transition zone. 

s 

The notations x' and y' sure used for the shock surface curvilineeu: 
orthogonal coordinates in Fig. 1. By using Stokes assumption and applying 
the hypersonic thin layer approximation, the governing equations for 
the shock transition zone can be e:q>ressed as [54] 

Continuity: 

p*V* = p*v* (3.13) 

09 a> S S 


X* -momentum: 

2 

p* + p*V*v* - C4/3)u* 3v*/3y* » p*V* + p* (3.14) 

^ oo 00 g g “* co as -^oo 

y ' -momentum • 

p*V*u - p* 3uV3y* = p*V*U* (3.15) 

00 00 ' 00 00 oo 

Energy : 

p*V*H* - Cjj*/Pr) 3/3y{H- (l-Pr)u*V2 - 

OO 00 

[l-(4/3)Pr]v*V2} = p;V*H* 
where H* * h* + (u*^ + v*^)/2. 

At the down stream edge of the trcinsition zone both v and (4/3) 

[p C3v/3y) ] are reduced to high-order quantities under a high shock com- 
pression ration. Therefore, a set of modified transport boundary condi- 
tions immediately behind the shock can be written as 

p*V* = p*v* 

00 OO ^ 


(3.17) 


(3.18) 


p*V*u* -y*(3uV3y) 
* * s s . 


P*V*U* 


P*+P*V*V* =a p*+ p*V*^ 

■*' » ^ 00 Q> 00 


p*V* (H*-H*)={ (Vi /Pr) [H*- d-Pr)u*2/2] } 

00 00 g CO g ■ * g ■ 


(3.19) 

(3.20) 


By introducing v » v sin a, u « v cos a and nondimens ionalizing 
all the quantities, the final modified Pan)cine<-Hugoniot conditions are 
obtained as 


p V* 

s s 


sin a 


u» = cos a - (e^ y /sin a) Ou'/Sy') 
s s s 

p ~ p + sin^ a (1-1/p ) 
s “ s 

h - h - (e^ y Pr sin a) Oh/Dy') + 
S 00 s 


(3.21) 

(3.22) 

(3.23) 


(1/2) [u' - cos a)^^ + sin^ a - v'^ 
s s 


u 


u' sin (a + 3) + v* cos (a + 3) 
s s 


V = V' sin (a + 3) - v* cos (a + 3) 
s s s 


(3.24) 

(3.25) 

(3.26) 

As mentioned before, use of slip boundaucy conditions are made in . 
investigating the shock layer flow phencniena at relatively high entry 
altitudes. 
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4. THERMODYNAMIC AND TRANSPORT PROPERTIES 
ThttxnodynaBiic properties for specific heat, enthalpy, and free 
energy, and transport properties for viscosity and thermal conductivity 
are required for each species considered In different flow regimes. For 
the precursor zone as well as shock layer, the general esqpresslon for 
total enthalpy, specific enthalpy, and specific heat at constant pres- 


sure are given respectively by 
• h + (u^ + v^)/2 
h « iC^h^ 

S * ^^iSi 


(4.1) 

(4.2) 

(4.3) 


However, specific relations for H and Cp for the tvio regions are quite 
different. 

For the precursor region, the relation for the specific enthalpy Is 
obtained by following the procedure described by Smith [1] as 

(4.4) 

where D and l represent the dissociation and Ionization energy respec- 
tively, and their values are available in [55] . The derivation of 
Eq. C4.4) essentially follows from the consideration of Eq. (4.2). If 
It is assumed that the Internal energy of each particle can be described 
only by translational and rotational inodes, then the relation for specific 


h" - 1.4575RT + (0.75RT + D)C + (1.25RT + 1/2)C 

H Ht 


enthalpy of each species can be expressed as 
ha. - §• RT + p/p - I RT 
hn - I* RT + I RT + p/p - y RT 

hRj - I* RT + I RT + p/p + I - j RT +I 


(4.5) 

C4.6) 

(4.7) 
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(4.8) 


3 5 

h„ = — RT + p/ + D * RT + D 

n 2 A 

h * I- RT (4.9) 

e 2 

Also, from the conservation of charged particles one cam write 


C /M 
e e 




(4.10) 


Now, for 85 percent and 15 percent He on volume basis (or 76 percent 
cind 24 percent He on mass basis), Eq. (4.2) is written as 

(Cc/Ma)ho^ * (0.26/4) C5RT/2) + ((0.74 - - C^)/2] (7RT/2) + 


[(5RT/2 +D))C„+ (7RT/2 + I) (C„ /2) + (5RT/2) (C /2) (4.11) 

H HJ HJ 


A simplification of the above equation results in Eq. (4.4). 

In the shock layer region, Eqs. (4.2) and (4.3) axe used to 
calculate H and C^. With representing the mole fraction of the ith 
species, the expressions for h^ emd are found from Refs. 56 and 57 
as 


H^ = RT[a^ + (a.^/2)T + (a^/3)T^ + (a^/4)T^ + 

(a5/5)T^ - a^/Tl (4.12) 

= R(a^ + a^T + a^T^ + a^T^ + a^T^) (4,13) 

where R is the universal gas constant C*l. 98726 cal/^ole ” K) and T 
is the local fluid temperature in K. For different species, values of 
the constants a^, a^# . . . a^ are given in [57], and for species under 
present investigation they are listed in Tal>le 3. It should be pointed 
out here that in this study, instead of employing Eq. (2.36b), Eqs. (4.1), 
(4.2), and (4.12) are used to calculate the enthalpy variation in the 
shock layer. This is because slightly better results are obtained by 
using the above set of equati<xis. 
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Table 3 Coefficient for evaluation of the specific heat at consteint 
pressure and enthalpy for various hydrogen/helium species. 
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For the shock^layer gas, the mixture viscosity and thermal con- 
ductivity are obtained by using the semi-empirical formulas of Wilke [58] 
as 


N N 

u “ I [x.M./iZ X.(|>. .)] (4.14) 

i-1 ^ ^ j-1 3 

N N 

K*Z ix.K./a x.4>..)] (4.15) 

1=1 j=l •* 

where 

♦ - [1 + Cli./U.)^'^^ CM./M 

* J J * J 


and is the moleculatr weight of species i. For hydrogen/helium 
species, specific relations for viscosity and thermal conductivity aore 
given in [59,60] . The viscosity of and He, as a function of tempera- 
ture, can be obtained from reference [59] as 

= (0.66 X 10“^) (T)^'^V(T + 70.5), N sec/m^ (4.16) 

“2 

u = (1.55 X 10 “®) (T)^'^V (T + 97.8), N sec/m^ (4.17) 

n0 

The thermal conductivity of and H are obtained frcm Ref. 60 as 

Kjj = 3.212 X 10 “^ + (5.344 x lo“^)T (4.18) 

K = 2.496 X 10 “^ + (5.129 x 10"®)T (4.19) 

H 

The viscosity of H amd thermal conductivity of He aore obtained from the 
relation between viscosity and thermal conductivity of monatomic gases 
as given in Ref. 58 by 

K = (15/4) (R/M)u (4.20) 

Very little information is available on transport properties of other 
species such as H^, e , etc. Fortunately, transport properties are 
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Important only in the boundary-layer region where the temperature is not 
high enough to produce these species. 

It should be noted that all relations presented in this section 
are expressed in dimensioneLL form. 

The heat transfer to the wall due to conduction euid diffusion is 
referred here as the convective heat flux euid is given by the relation 
CIS, 46] as 

/ ” 

q « -e2[K(3T/9y) + (w Le/Pr) Z OC,/3y)h,] (4,21) 

where Pr is the Prandtl number, Le is Lewis number and the value for Le 
is taken to be 1.1 [45,46] in the present study. 
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5 . CHEMICAL REACTIONS 


Analyses of chemically reacting flows are usually simplified by 
assuming the chemical equilibrium behavior of the gas mixture. While 
this assumption may be justified in some cases, in many realistic prob- 
leaa this may lead to serious errors. Thus, in order to understemd the 
degree of physical reality, it becomes essential to analyze the complex 
gas mixture under the conditions of chemical nonequilibriinn . In this 
chapter, information on chemical equilibrium and nonequilibrium reactions 
and reaction rates are provided for the shock layer gas mixtxire of a 
Jovian entry body. 

5.1 Chemical Equilibrium 

In the chemical equilibrium case, a computer code developed by 
Sutton [26] is used in this study. The niamber density of eight chemical 

^ • "I" ^ 

species, H^, H, H , H , e“, He, He and He are calculated by the chemi- 
cal reactions euid rate constemts given in Table 4. In general, the 
reactions can be expressed by 

L a. A, -t- Z b.B. (5.1) 

X i -► XX 

3 

The number density of particles (particle/m ) is related to the 
equilibrium rate constant and can be expressed as [28] 

K. - [nN^i(B.)l/[nN^i(A.)] (5.2) 

3 X X 

The conversion equations for hydrogen and helium nuclei and 



(5.3) 

(5.4) 



(5.5) 
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Table 4 Reaction scheme and rate constants for 
chemical equilibriiam conditions. 


1. 

2 . ♦ •* 

3 . + •” 

4. »^^*am** + •" 

5. ♦ •" 


BATE COWSTMITS PARtlCLES/M^ 

ln^ • 4.699C22 <l-«xp(-6331/T) exp(-51964/T) 


- 2.411E15 


•xp (-157810/T) 

- 9.64SE1S 


•xp(-2aS287/T) 

- 2.411£15 

T 

•xp (-631310/T) 

- 9.643E1S 

T 

mxp (-87S0/T) 


Table 5 Reaction scheme and rate constaints for 
chemical nonequilibrium conditions. 



AMCtiaiia 

Bata constants in cm^ sec ^ i 

I. 

H + • j H* a« 

1/2 

kj^-2.27E13 exp(-LS78E5/T^) 

3. 

a« • • ^ + 2a 

lt2“l*33El3 axp(-2.8S2E5/T^) 

3. 

H ♦ a j •* 

«• ♦ • 5 * + 2a 

kj- 4.11E13 axp(-L.160ES/T^) 

4. 

■a 4’ a ; Ha* ♦ a, 
■a* ♦ a J Ha* • 2a 

k^- 2.24E13 axp(-2.320E5/T^) 

S. 

« 4 « J ■* 4 H, 
B*4H*H 4«4H 

kj-6.20E10 ajq>(-1.160E5/T) 

«. 

H 4 Ha 4 H* 4 Ma, 
R*4Ha4H 4a4ila 

kg-4.89E10 a^(-1.160E5/T) 

7. 

H2 4Ba4H4H4Ha 

k^-4.33E18 t”^ [l-exp(-15E8/T^) I 



axp(-52340/T) 

• . 

»a ♦ ; H 4 H 4 

k «2.S k, 
8 7 

9. 

4H4H4H+H 

kg-14.0 k^ 

10. 

4 4 H 4 H 4 

"lo’S 

11. 

Hj4a4H4H4a 

‘u"'9 
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The number densities of the hydrogen and helium nuclei are calcu^ 


lated by 




(5.6) 

C5.7) 


where 

M - 2.016x + 4.003X 

o He 

In the above equations, represents Avogadro's const£uit^ p Is the 

3 

mixture density in g/cm , x is the mole fraction of moleculeir hydrogen^ 

arid is the mole fraction of helium. 

He 

The solution procedure for obtaining the eight \inknown number 
densities is discussed in [26] . The closed-form solutions are obtained 
by solving Eg. (5.2) for each reaction independently. This is accom- 
plished by setting the appropriate values in Eqs. (5.3)- (5.5) equal to 
zero if the species are not present in the reaction. The closed-form 

3 

solutions for the number densities (in particles/cm ) of each species 
are given by 

“ 2 ' ^H ^ (K^/8)U + 8N^K^)^/^ - 1] 

h'^; N^+ = ^^/2)U1 + - 1] 

H: 

<= = ®l/^l [(1 + . ij, D 3 - K 3 + 

= CD3/2) [(1 + - 11 , D3 - + Ng + 


He: N = N - N„ + N ++ 

He He He He 


e"; N_=N++N++2N++ 

e H He He 


H-: V = Ve-/^5 


C5.8) 
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5.2 Chemical Nonequilibrium 

When chenical reactions proceed at a finite rate, knowledge of the 
rate of production terms, which appear in the species continuity 

equations, are required. The reaction scheme describing important 
colli sioned. processes in hydrogen-helium ionizing shock waves has been 
modeled by Leibowitz [8] after the results of argon ionization studies. 
Eleven separate reaction steps describe the dissociation of molecular 
hydrogen, excitation of electronic states of hydrogen and helium, and 
ionization of the atcmic hydrogen and helium by collisions with atoms 
and electrons. A complete discussion on these reaction schemes is avail- 
able in [61] euid [62] . The eleven reactions and corresponding rate con- 
stants aLtre given in Table 5. 

In a ccmiplex gas mixture containing a total of i species, of which 
are capable of undergoing m elementary chemical reactions, the chemi- 
cal equation for the general elementary reaction r can be written as 
[28] 


Z Z 

Z a. X. Z b. 

r,r X •*- . 


i=l 


i=l 


X. 
i,r 1 


(5.9) 


where a, and b. are the stochiometric coefficients appeeuring on 
1 ,r f 

the left and right in the reaction r. By applying the principle of 
detailed balancing, the backward rate constant, is obtained by 

dividing the forward rate constant, K_ , by the equilibrixmi constant 
K which is given in Table 1, 

C /I" 

The total rate of change in X^ is given by the relation [28] 


dx.' m Z ^ 

3~= Z (b. -a. )K. ^tn (X,)®l,r - i ",(X.)‘’i,r] 

dt x,r i,r f,r i K i=l x 


(5.10) 


i=l 


c,r 
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This is the general rate equation for a complex gas mixture. The rate 
of production of chemical species, now can be expressed by 

= M*(dxVdt) (R*/P*V*) (5.11) 

X X i n * 

Equations (S.'S)*- (.5.11) , along with other fluid mechanical equations, 
equation of state, and the electron energy equation, are solved numeric** 
ally to obtain the concentration of all species. In order to have a 
reasonable rate of convergence in the numerical scheme, however, it is 
important to express the rate of production term in a proper form. This 
is accomplished by splitting into two sepeurate contributions as 


[15,63,64] 

^./p = (ft.)® - (ft.)^C. (5.12) 

X X XX 

The reasons for doing this are explained in the cited references. 

5.2.1 Electron Temperature 

Because of a large ratio of atom (or ion) mass to electron mass, 
electrons transfer energy rapidly by collisions with other electrons but 
only slowly by elastic collisions with atoms or ions. Conseq[uently, a 
different temperature is given to atoms (heavy particles) and electrons 
in the same gas. The electron temperature is obtained from the solution 
of the electron energy equation. A detailed discussion of the electron 
energy equation is given by Appleton and Bray [65]. For a one-dimensional 
steady shock wave in a H^-He mixture, the resulting equation can be 


expressed as [61] 

3 [el«^VR(T-T^)-0^(R^-R^^+R3-R3^)-e^^ <*‘2-'‘2r*V^4r >=° 


(5.13) 


where 


^ek = \ ^e 2ek 


(5.14a) 
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(5.14b) 


V - (8 k T /Tim ) 
e e e 


1/2 


Qgjj * Aj^[exp(-aj^E)-exp(-bj^E)] 

E * (1/2) 


V ■ 4(kT /2 ttv) 
e 


1/2 


U » MaMb/(Ma+Mb) 


* I '^ek'^k 
k 


(5.14c) 

(5.14d) 

(5.14e) 

(5.14f) 

(5.14g) 


In Eq. (5.13), le] represents the ccaicentration of electrons, 8 ^^ is 
the ionization energy per mole of species k, smd and eu:e the 
forward auid backward production rates for electrons respectively. In 
Eq. (5.14b), represents the collisional frecpiency, n^^ is the number 
density of species K, is the average electron velocity, and 
represents the elastic collision cross section for species k. In Eq. 
(5.14e), E represents the relative kinetic energy, V is the relative 
speed, amd \x is the reduced mass. The values of coefficients a^^, 
and b appearing in Eq. (5.14d) are available in Ref. 59. By substituting 
Eqs. (5.10) and (5.14) into Eq. (5.13), an explicity expression for the 
electron temperature, in terms of the species concentration and heavy 
particle temperature, can be obtained as 

- t(Vk3)(X2-X^/K^^„) + 0C2+k4)(X3-X5/K^^^^)]/X^ 


(5.15) 


wnere k. , k* , k. , k . eure rate constamts in Table 5 and 
1 2 3 4 

X, = l/(3[e]m vR) 

1 e 

*3 - 

X 4 = ejj[H+] t'e] 


’‘5 


(5.16a) 

(5.16b) 

(5.16c) 

(5.16d) 

r5.16e) 
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Use of the electron temperature is made in evaluating the radiative 
flux in the shock layer. 
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6. RADIATION MODEL 


An appropriate expression for the radiative flux, q , is needed for 
the solution of the energy equation presented in the second chapter. 

This requires a suited^le tremsport model emd a meaningful spectral model 
for variation of the eUasoxption coefficient of the gas. In this chapter, 
appropriate es^ressions for the spectral and total radiative flux are 
given and information on the spectral aibsorption by the hydrogen/helium 
gas is presented. 

6.1 Radiative Flux Equation 

The equations for radiative transport, in general, are integral 
equations which involve integration over both frequency spectrum and 
physical coordinates. In many physically realistic problems, the com- 
plexity of the three-dimensional radiative transfer can be reduced by 
Introduction of the '*t£Uigent sleds approximation.** This approximation 
treats the gas layer as a one-dimensional slab in calculation of the 
radiative transport. Radiation in directions other than normal to either 
the body or shock is neglected. Discussions on the validity of this 
approximation for planetary entry conditions are given in [66-70] . 

As mentioned eeurlier, the tsuigent slab approximation for radiative 


transfer is used in this study. It should be pointed out here that the 
tangent slab approximation is used only for the radiative transport and 
not for other flow veuriables. For a nonscattering medium and diffuse 
nonreflecting bounding surfaces, a one-dimensional expression for the 
spectral radiative flux is given by [27-29] 
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where 


T = /^ a (y* )dy* 

V V ' 

o 

1 n-2 

E (t) = / exp(-t/y*)vi dy 
n 

o 

== (hv^/c^) [exp (hv/KT) -1] 

The quantities B^(o) and represent the radloslties of the body 

surface and shock respectively. 

The expression of total radiative flux is given by 

<30 

q„ = / q„ )dv ( 6 , 2 ) 

RV V 

0 

To obtain specific relations for the total radiative flux for the pre* 
cursor and shock-layer regions, it is essential to know the spectrel 
absorption characteristics of the absorbing-emitting species in these 
regions. 

In the precursor region, the radiative contribution from the free 

stream usually is neglected. For a diffuse, nonreflecting shock front, 

the expression for one-dimensional radiative flux for this region is 

obtained from Eqs. (6.1) and (6.2) as 

q^(n) = 2 /“ (q^(O)E^ (K^n) + 

0 

TTK /“ B (T)E^[< (n-n ’ ) ] dii • }dv (6,3) 

V V 2 V 

0 

where q (0) = G^irB^(T^). In obtaining the above equation, it was 
assumed that the absorption coefficient ic^ is independent of position. 
The information on the spectral absorption model for hydrogen/ 
helium species in the precursor region is given in [42] and is briefly 
discussed in subsection 6.2, The model essentially consists of 
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approximating the actual absorption of active species by three differ- 


ent step models. For this model, Jq. (6.3) can be expressed as [43] 
q (n) = 2 tt E { (15/7T^)q(0)E (k .n) f ^ [v^/(e^ - l)]dv 


i«l 


n 


3 i 


'’21 


’'ll 


+ K. f E 2 [K.(n - n')l / B (T)d dn'} 

'’li 


(6.4) 


where V = hv/kT^ amd q(0) = eaT\ In writing the above equation, it has 
been assigned that the shock front radiates in the precursor zone as a 
gray body . 

In the shock layer, the radiative energy from the bow shock usually 
is neglected in comparison to the energy absorbed and emitted by the gas 
layer. This implies that the transparent shock front does not absorb but 
emits radiation. The expression for the net radiative flux in the shock 
layer, therefore, is given by 


q„ = 2 / [q (0)E, (T ) + f B (t)E_(r - t)dt 

R V 3 V V 2 V 

0 0 


/ B (t)E^ (t - T )dt]dv 
V 2 V 

T 

V 


(6.5a) 


In this equation, the first two terms on the right represent the radia- 
tive energy trcinsfer towards the bow shock while the third term repre- 
sents the energy transfer towards the body. Upon denoting these contribu- 

4* * 

tions by q and q , Eq. (6.5a) CcUi be written as 

R R 




(6.5b) 


A few spectral models for absorption by the hydrogeh/helium species 
in the shock layer have been proposed in the literature [22-26] . For, 
Jovian entry conditions, the absorption by helium is usually neglected. 
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The spectral absorption of hydrogen species was represented by a 58-step 
model by Sutton [26] auid was approximated by a 30-step model by Tiwari 
euid Subramanian [27] . The results of these step models aure ccoipared with 
the detailed model of Nicolet [22] in [27]. The 58-step model proposed 
by Sutton is employed in this study. The details of radiative edssorption 
and computational procedxire aure given in [26] . The information on spec- 
tral absorption by this model is summarized in subsection 6.2. In 
essence, the step model replaces the frequency integration in Eq. (6.5) 
by a summation over 58 different frequency intervals. In each interval, 
the 'absorption coefficient is taken to be independent of frequency. For 
this model, Eq. (6.5) cein be expressed as 
N 

q„ = 2ir E {e B (T )E. [ (y')d '] 

.-vvw3 v”" y 

3=1 0 

y y 

+ / a^(C)B^(4)E2C / a^(y')dy*]d€ 
o 5 

- r a^(C)B^(C)E2[ f a^(y')dy‘]dC} (6.6) 

y ^ ^ y ^ 


where y^ denotes the shock location 6md N represents the number of 
spectral intervals. In each of the jth intervals, the absorption co- 
efficient is assumed constant while the Plemck function is not. In 
accordcince with Eq. (6.5b), Eq. (6.6) Ceui be expressed in terms of q^ 

auid q and for a gray body one finds 
R 


+ N 

%(y) = (47Th/c2) E {eF(v.,TjE^[^/^ a^. (y')dy’] 
j=l 3 w o V3 

y y 

+ / (KT/h)'+ F(v ,T)a^. (^)E 2 [ / a (y')dy*Jd5} 

o J J ^ J 


(6.7a) 
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N Y ? 

q”(y) » - (4TTh/c^)2 { /^CKT/h)** F (v . ,T) .E^ I / a .(y*) 

j»l y 3 y VD 


X dy'ld?} 


where 

F(v.,T ) 
3 w 

F(v^,T) 


= /^^{vV[exp(hv/KT ) - l]}dv 

w 

''31 

=* /*' {v^/[exp(v) - 1] }dv, V » hv/KT 


(6.7b) 


From the knowledge of the temperature distribution normal to the 
body, Eqs. (6.7) CeUi be solved by numerical integration over frequency 
and space. The final temperature profile is obtained through an itera- 
tive procedure. Use of Eqs. (6.7) is made in obtaining the radiative 
flux towards the body and shock as well as the net radiative flux. 

For evaluation of the radiative flux, usually it is essential to 
express the exponential integrals simpler approximate forms. 

Quite often, these integrals are approximated by appropriate exponential 
functions [28,29]. In this study, it was established that better results 
are obtained if the exponential integrals are expressed in series form 
for small and large arguments. The series expansion of the exponential 
integral of first order is given as 
For t < 1: 

2 3 

E^(t) = - 0.5772 - £n t + t - I(3)T + • * * (6.8a) 

For t > 1: 


E^(t) = exp(-t) 


0 It 2t 


+ a 


3t + t‘ 


+ bi^ + b 


2t 


+ b 


3t 


+ t**) 


(6 -8b) 
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where 


* 0.26777343 » 3.958469228 

0 0 

=* 8.63476089 • 21.09965309 

a^ » 18.059016f73 b^ » 25.63295614 

a^ » 8.57322874 b^ - 9.5733223454 

Relations for exponential integrals of higher order are obtained by 
employing the recursion relations given dLn [29] . 

6.2 Radiation Absorption Model 

Appropriate spectral models for gaseous absorption are needed for 
solutions of the radiative flux equations. Information on spectral 
absorption by the precursor emd shock-layer species is presented in this 
section . 

6.2.1 Spectral Absorption Model for Precursor Region 

In the precursor region, the photoionization absorption coefficient 
is a continuous nonzero fxinction of photon energy (because of bound-free 
transition) for all values of photon energy that exceed the ionization 
potential of the atom. Similar remarks apply to the photodissociation 
and radiative recombination. A critical review of ultraviolet photo- 
cdDsorption cross sections for molecules of astrophysical emd aeronomic 
interest, available in the literature up to 1971, is given by Hudson [20]. 
Specific information on photoionization euid absorption coefficient of 
molecular hydrogen is available in [20,21, 71-74]. 

Photoionization and absorption cross sections of H^, as obtained 
from Refs. 20, 21, ^uld 71-74, cure plotted in Fig. 2. From this figure, 

O 

it is evident that the ionization continuiom starts at about 804 A auid 

O 

continues towaurds lower wavelengths. Between the wavelengths of 600 A 
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Fig. 2 Absorption cross section of H in ultraviolet region. 





and 804 A, the absorption cross section for the ionization continuian Is 
included in the total absorption (i.e., absorption due to ionization as 

o 

well as dissociation). For wavelengths below 600 A, however, the ioniza- 
tion continuum edssorption is equal to the total eOssorption. Hie total 

o 

lObsorption cross section for the continuum rzmge below 804 A can be 
closely approximated by the two rectangles (I auid II) shown in the 
figure with broken lines. The ratio of the ionization cross section to 
the total absorption cross section (i.e., the value of Y^) is taken to be 
unity for rectangle I and 0.875 for recteuigle II. For wavelengths greater 

O 

than 804 A (where h is below ionization energy) , the value of is taken 
to be zero. Little information is available in the literature on the 
absorption cross section for dissociation of molecules. Hiere is 

o 

strong evidence, however, that photodissociation starts at about 2600 A 

O 

auid continues towards lower wavelengths to 5d30ut 750 A [69,71]. Hiere 
au:e also a few diffuse bands in this spectral range [71,73]. Thus, it 
becomes difficult to evaluate the absorption cross section in this spec- 
tral range. For this study, the absorption cross section in the spectral 

O O 

range between 804 A and 2600 A was approximated by rectangle III. The 
specific values of ct(v) for the three rectauigles are found to be 
a^(v) =* 4.1 E-18, =8.2 E-18, and =2.1 E-18. The value 

of Y^ is taken to be zero for rectangle I auid 0.125 for rectangle II. 

The numerical procedure for employing this model in the radiative flux 
equations is discussed in detail in Ref. 42 aind is summarized in Chap. 7. 

6,2.2 Spectral Absorption Model for Shock-Layer Region 

As mentioned earlier, the 58-step model proposed by Sutton [26] for 
spectral absorption by the hydrogen species in the shock layer is employed 
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In this study. For atomic hydrogen, all three transitions (bound^bound , 
bound~free, and free-free) are considered. The total absorption of the 
jth step Is a summation of the average absorption for the 1th transi- 
tions In the jth step, l.e.. 


•'j " I "ij 

V +AV 

tCj^. ■ (1/Av ) /■* ^ (C^dv 


(6.9a) 

(6.9b) 


•c^ * f (T,N^,v) 

where represents the number density in cm 


(6.9c) 


For the free-free transition, the absorption coefficient is cal- 


culated by 
H 


tc„ - (2.61E - 35 )N N +/(v3t^/^) 
fr e H 


( 6 . 10 ) 

The adssorptlon coefficient for bound-free trsurisitions is calculated 


by employing two separate relations as 
H ^ 

%f “ - 14)(Nj^v3) I (l/n|)exp(C^), 1 1 ^ 4 (6.11a) 


Sf - <® ■ 31E - 20) (TN^v 3 ) exp (C^ ) ) , 5 < n^< n j ^ 

where 

- (-157780/T) 11 - (1/nJ)] 

Cj - (-157780/T) (1 - 6/13/6) 

- [ (157780/T) (1/25 - 6/13/6)] - 1 

6 « (1.79E - 5) (N^'^^)/(T^/^) 

e 

In the above ecjuations, n^ represents the principal quantum numbers, 6 
is the reduction in ionization potential in eV, and the values 157780 
and 13.6 are the ionization potential in K and eV respectively. 
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The bound-bound transitions are Included for principal quantum 
numbers up to five. The absorption coefficient is calculated by using 
the relation 

k” = SL(v) (6.12) 

on 

where S is the line strength and L(v) is the line shape factor. The 
line strength is given by the relation 

S = (l.lOE - 16)fn2N„ expl (-157780/T) (1 - l/nj)] (6.13) 

The line shape factor is given by the relation 

L(v) = Y/{TT[y2 + (v - (6.14) 

where is the frequency at the line center and y is the line half- 
width, and these are given by 

V = 13.6[(l/n2) - (l/n^)] (6.15) 

0 X. U 

y = a[1.05E 15 (n^ - n2)N^/^] (6.16) 

u X. e 

The constant a in the above equation is taken to be 0.642 for the first 

line and unity for the remaining lines. 

The absorption coefficients for the free-free and bound-free transi 

tions of the negative hydrogen are 

= (6.02E - 39)N N /v^ (6.17) 

f r n e 

= (2.89E - 17) (e'* - 4^3 + 3.64p2 + 0.73B)N - (6.18) 

where B = 1.502/v. The threshold for the bound-free transition of H"" 
is 0.757 eV. 

The absorption coefficient for molecular hydrogen in the jth step 
is obtained in accordance with Eq. (6.9) cind is expressed as 

ic^2 == f . (T)N _ (6.19) 

where f^ (T) is dependent on the particular step. The molecular bands 
cover the steps from 7 to 17 eV. 
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Further details on constructing the step<<^f unction model emd utilize 
ing it in the radiative flux equations are given in Refs. 25'-27, 
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7 . SOLUTION PROCEDURE 


An iterative procedure has been used to couple the precursor and 
shock-layer solutions. In this method, the shock-layer solutions are 
obtained first with no consideration of precursor effect. From this 
solution, the radiative flux at the shock front (which influences the 
precursor region flow) is determined. By employing this value of the 
radiative flux, different precursor region variables are calculated 
through use of Eqs. (2.18) through (2.25). Values of these flow variables 
are obtained just cihead of the bow shock, and then the Rankine-Hugoniot 
relations are used to determine the conditions behind the shock. These 
conditions are used to obtain new shock-layer solutions from which a new 
value of the radiative flux at the shock is calculated. The procedure 
is continued until the radiative flux at the shock becomes invariant. 

The solution procedures for the precursor and shock-layer regions 
are described in some detail in following subsections. 

7.1 Precursor Region Solution 

As pointed out earlier, two methods (the small perturbation theoiry 
and the thin-layer approximation) aire employed in this study to investi- 
gate the precursor region flow phenomena. The solution procedure for 
these methods is discussed separately in this section. 

7.1.1 Small Perturbation Theory 

Since the problem treated by thin layer method is linear, it is 
permissible to obtain a solution for arbitrary frequency, and then 
integrate this solution over the spectrum to obtain the general solu- 
tion. Thus, in the development that follows, flow-field perturbations 
will be considered for a unit frequency interval. Consequently, Eqs. 

(2.11) and (2.12) now can be written as 
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C7.1) 


3H_ /3z » N a(v)/(p v3) 

T H oo CO \) 

1 “2 


3C /9z * fm, aCv)/(p V hv)]H 

n^+ 1 H2 I to CO V 

3C y92 “ Cm a(v)/(p V hv)]H 

IT 1 D '^00 00 V 


C7.2a) 

C7,2b) 


where Y^^ and Y^ represent photodissociation amd photoionization yields 
respectively. 

It cam be shown that the flow under consideration is irrotational 

[1,2]. Thus, there exists a potential (}) such that 

(7.3) 

For z-direction, integration of Eq. (2.8) results in 

(7.4) 


= n 


“ -Cym2) 34>/3z = -(yM^)V^^ 


Eq. (2.9) ccin now be expressed as 
+ 3p^/3z = 0 


(7.5) 


In order to evaluate 3p^/3z and to relate to other variables, it is 
necessary to consider the gas model and radiation. For the precursor 
region gas model, the expressions for pressure and enthalpy variations 
are given by Eqs, (2.14) and (2.15) respectively. Now, in order to ex- 
press the governing equations in terms of perturbation potential, first 
p^ is eliminated by combining Eqs. (2.14) cind (7.4), The resulting 
equation is then differentiated with respect to z and use is made of 
Eq. (7.2). Next, Eqs. (2.15), (2,16), (7,1), and (7.3) are combined to 
give 

3p,/3z = -r 32(|)/3 z2 -PH (7.6) 

1 V V 

where 

r = 0.727 Y m 2 (7.7a) 


P = a + b /hv 

V V \) 


(7.7b) 
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a - N a(v)/(p„ H J (7.8a) 

V H2 * “ * 

b * -(a m,/2)I(I - 0.89 RT )Y^ + (2D - 1.89 RT ) Y 1 (7.8b) 

V V 1 « I * D 

Upon combining Egs. (7.5) auid (7.6) , the governing equation for the flovr 
is obtained as 

7 ^ < 1 > - r B^6/dz^ -PH (7.9) 

x,y^ ^ V V 

For the eocisymmetric case, this is e^q>ressed as 

r“^ 3/3r (rd<t>/dr) - r B^<f>/dz^ - (7,10) 

Eqs. (7.9) €uid (7.10) 2 u:e seen to be the classical potential equa- 
tions for compressible flovr with a forcing term proportional to radiation 
added. The potential for the flow induced by a radiant source with a 
spectral distribution is obtained by integrating the contributions of 
each frequency as 

* = (|> dv (7.11) 

0 

As discussed by Smith [1,2], solutions of the governing equations, 
presented in the previous section, C 2 ui be obtained in special cases 
depending on the model used for the distribution of spatial radiation. 

If the radius of the radiating gas cap, R^, is large compared to the 

photon mean free path, then the problem can be treated like radiation 

from a plane source. On the other hauid, when the radius of the radia- 
ting gas cap is small, then the problem caui be treated like a spherical 

point source for radiation from the gas cap and a cylindrical point 

source for radiation from the wake. Note that, in general, R^ may not 

be the same as the radius of the bow shock, R . 

s 

7. 1.1.1 Radiation From a Plane Source . For radiation from a plane 
source, it is essential to integrate the contribution over the plauie. 
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as attenuated by passage through the absorbing medium. The relation for 
in this case, is given by [28] 

C7.12) 

where < 1 ^( 0 ) is the spectral radiative flux density at the shock, wave, 
is the spectral aUbsorption coefficient, auid is the exponential 

integral of order n . The expression for (which may also be interpreted 
as inverse of the photon meaui free path) is given by 

K « N <T(U) C7.13) 

“2 

In this form represents the absorption coefficient of molecules. 

If the number density (and hence, k^) can be taken to be independent 

of z (which is a good approximation for small ionization and dissocia- 
tion) , then the optical depth is defined’ by 

C “ (7.14) 

For the plane radiating source (where 7^ * 0) , therefore , a combina- 

tion of Eqs. (7.9), (7.12), and (7.14) results in a simpler expression, 
the integration of which results in [42] 

((> - -[2 q^(0)/(rtc2)]E^(-C) (7.15) 

where the boundary condition of (3(j5/3?)-^0 as has been used. 

Prom Eq. (7.3), the velocity perturbations, ahead of the shock 

front, can now be written as V_ ** V_ * 0, and 

lx ly 

V, = -[2 P, q (0)/(rK:^ )]E,(-?) (7,16) 

iz V V V 3 

From Eq. (7.4) , the expression for pressure perturbation is foTind to be 
Pi “ [2 Y q^(0)A^]E3(-C) (7.17) 

where it was assumed that Similaurly, the expressions for 

density perturbation, total enthalpy, static enthalpy, and species 
concentration are found to be [2,42] 


57 



radiates as a cylinder of infinitesimal radius emd length. For both 
cases, the incident radiation at 2 uiy field point s is given by [1,2] 

* (A^/s^) exp(-K^s) (sin 6)^ (7.22) 

In this equation, represents the radiative strength of the source, 
s is the distamce from the source euid 6 is the emgle between the free 
stream velocity vector and a line from the field point to the center of 
the source. The superscript j is equal to zero for a spherical point 


source euid one for a cylindrical point source. 

Eq. (7.22) c£ui be substituted in Eqs. (7. 9) •(7.10) to obtain the 
corresponding equations for the perturbation potential. Within the con> 
fines of the assumptions made in obtaining Eq. (7.22), however, both 
problems (spherical as well as cylindrical point source) can be con* 
sidered to be axisymmetric . The governing equation for the perturbation 
potential, therefore, can be written as 
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Ui 
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Fig. 3 i^proxination of radiating flow by point sources. 



C7.23) 


n”^ a/3n(Ti3(|) V8n) - r =• a (sin 9)^ 

where 

u = KS, n = Kr, A*A(v1=k P A 

V V V vvv 

A procediire for '“general solution of this equation is suggested by- 
Smith tl] . For entry flow, however, M^>>1 and Eq, (7.23J can be solved 

OP 

by e:q>euiding <j> in a series in (l/F) in the vicinity of the body. Thus, 
one can express as 

<J)’ - -(Kyr)[F^(c,n) + (i/r).F^^^^a,Ti) + d/rj 

+ . . .] (7.24) 

where F^'s are function for perturbation, potential. Substitution of 
this relation into Eq. (7.23) gives 

3^F^/3C^ = U ^ e:q>(“li) (sin 6)^ (7.25) 

cuid 

^2p (n)/3^2 , 3/3n(Ti 3 f/^“^V3ti) (7.26) 

J 

The problem, therefore, is reduced to quadratures in the vicinity of the 
body. In the present analysis, only the terms in (l/F) will be retained. 
By integrating Eq. (7.25) twice, the e:q>ression for F^ is obtained as 

F. (?,n) “ exp(-y ) (n/y )^ (C • ? )d? (7,27) 

i o 0 0 0 0 

^ -»00 

2 

where y^ ~ For convenience, let us denote 

G. (C,T 1 ) = 3F./3C » y'^ exp (-y^) (n/y ^)^ d;^ (7.28a) 

H. (C,n) = 3F./3n * y ^ e:q)(-y ){(n/y)^ [n + 

3 3 ® 

(n/y^) (2 + j)l - j }(? - (7,28b) 

with these definitions of F ^ , G ^ , and , the perturbation quantities 

can be expressed as [2,42] 
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(7.29) 


♦ * « -(A/r)p^(i;,Ti) 

^Xr “ (7.30) 

Vi^ - -(A/DG^ (C,n) (7.31) 

p^*YAGj(C,n) (7.32) 

- (A/DG^(C,n) (7.33) 

- (K^ \/p^ v3)G^(C,ti) (7.34) 

Cjj - (mj^ k2/p^ hv)Yj^(v) G^(C,Ti) (7.35a) 

^ 2 + " \ ^l/Pco hv)Y^(v) Gj(J;,n) (7.35b) 

Note that for the case of spherically radiating point source j » 1 in 


the adx)ve equations. Also, these equations are obtained for arbitrary 

frequency. The expression for total potential, for this case, can be 

obtained by c<xnblning Eqs. (7.11) euid (7.24) as 

♦ - -Cl/D /“ [A(v)/icv]F (C,n)dv (7.36) 

0 ^ 

Furthermore, it should be noted that the above solutions are valid in 

-2 i 

the region where [y exp(-y) (sin 6)*^] does not vemish. This is the 
case of spherically symmetric flow 2 diead of the entry body and is of 
primary concern in the present study. Other cases Involving cylindrical 
point source are discussed in [1,2] . 

The procedure for expressing the perturbation eq[uations in terms of 
the photoadssorption model employed in the precursor region is given, in 
detail, in [42]. 

7.1.2 Thin Layer Approximation 

A direct integration of Eqs. (2.23) through (2.27) results in the 
following governing equations for the precursor region 
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(7.37) 


pv = p V 

^ OO OD 


p . V (u - u ) = 0 

CO 00 


P«o + CP - p«> 


P« v^ CH - + ‘Ir “ 0 

p« ° 

00 00 Q Ot 

where it has been assumed that * 0. 

R« 


(7.38) 

(7.39) 

(7.40) 

(7.41) 


In Eq. (7.40), H represents the total enthalpy and Is given by a 

combination of Eqs. (4.1) amd (4.4). The expression for the radiative 

flux, q , is given by Eq. (6.4), For the present application, Eq, (7,41) 
R 

will be written for atomic hydrogen and hydrogen ions. By following 
the procedure outlined in [1,42] the expressions for species concentra- 


tion are found to be 
N 


- 1 . 


C = 2B. Z E^(T,)(kT^) l(v.) 

H 4.,D. 3i S 1 

1=1 1 


N , 

»+ - 23. Z E, (T.XkT ) I(V?) 

^^2+ ^ i=l ^i ^ ^ ® ^ 


(7,42) 


(7.43) 


where 


3, = (IS/TT*^) [q(0)m,/(p V )] 

4^ ^ 00 00 

V 

I(v?) = {vV[exp(v) - l]}dv 

^li 

and m^ represents the weight of the molecule in grams per molecule. 

Note that there are nine algebraic equations to evaluate the nine 

unknowns p, v, u, p, h, H, C, C . The solutions of this set of 

H H 2 + 

algebraic equations are obtained by using the Gauss-Seidel method [74] . 
The properties at the infinity are used as the first initial guess in 
the Gauss-Seidel method. The iteration is continued until all the 
quantities in this region become invariant. The flow chart of the 
computational procedure is illustrated in Figs, 4 cind 5, 


62 














Pig. 5 Plow chaurt for subroutine PERC used in the precursor region 
solution procedure. 
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7.2 Shock- Layer Region 


A numerical procedure for solving the viscous shock-layer equations 
for stagnation and downstream regions is given by Davis [14] , Moss 
applied this method of solution to reacting multicomponent mixtures in 
[15,17], A modified form of this procedure is used din this study to 
obtain solutions of the viscous shock-layer equations. In this method, 
a transformation is applied to the viscous shock-layer equations in order 
to simplify the numerical computations. In this transformation most of 
the v^u:i^d^les eure normalized with their local shock values. The trauis- 
formed variadsles are [15] 


n - y/y^ 

p - P/P 3 


5 - X 

p" ' P/P3 

K = K/K 

s 

u » u/u 

T = T/T 


s 

s 

p p ps 

— » v/v 
V ' s 

H - H/H 

s 

(7.44) 


The transformations relating the differential quantities are 


auid 


^ {dy^/d5)n^ ( ) 


^s 3y2 y2 8n2 


1 


( ) 


(7.45) 


(7.46. 


The trauisformed equations cam be e}q)ressed in a general form as 

32w/3ti^ + a^3W/3n + a^W + *3 + a^3W/3C = 0 (7.47) 


The quamtity W represents u in the X-momentum equation, T in the tempera- 
ture energy equation, H in the enthalpy energy equation, and in the 
species continuity equations. In most cases, the coefficients a^ to a^ 
to be used in this study are exactly the same as given in [15]. However, 
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there is one exception. Since radiation Is included in the present 
study, the coefficients of the energy equaticm are different frcsi those 
used in [15]. For exai^>le, in the enthalpy energy equation, coeffic* 
ients » and a^ are the same as given in [15] , but a^ is different , 

and this is given by 


r,i 

a -— «■ 


P y^ 
3 r s 


y yH 
s s 


( £_ 

y_ 3 ti 1 + yr\< 

3 S 


cos 9 


r + y n cos 0 

s 




y P_p V V 
3 ars 3 

e^M yH 


. y P P , 3q„ 

3p 3 r r,s - 1 ^ 

3ti ~ 2 IT “ y 

e^y H y ■*s 


s s 


+ 


R 1 + YgHK 


COS 6 

r + y cos 0 


01 


where 


N 3C, 

ys . O.ly . j 

y Pr,s^ Pr . , ”i 3n 

3 1*1 


Pr 




T1 y2 Kyu^ 
S 3 

1 + n«: 


( 7 . 48 ) 


Other transformed equations are the same as given in [15] • 

The surface boundary conditions in terms of transformed variables 


are 

il=0, v = 0, T*T (7.49) 

w 

The trcois formed shock conditions are found to be 

u = v = T = ii»p*p*l (7.50) 

at n * 1. 

The second order partial differential equations as expressed by 
Eq. (7,47), along with the sTirface boundary and shock conditions, are 
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solved employing an Implicit flnlte<^lfference method. In order to 
^taln numerical solutions for the dovmstream region, It Is necessary 
to have an accurate stagnation streamline solution. Since the shock 
shape Is affected by the downstream flow, a truncated series of shock 
standoff distance Is used to develop the stagnation streamline equations. 
As such, the shock standoff distance Is eiq^ressed by 

y. “ ^1. + ^23 

since C Is small and the curvature k Is approximately, one in the 
stagnation region. It Is logical to say that (see Fig. 1) 

0 » C (7,52) 


Since 9 « (yi/2) -6, there is obtained 

a « 9 + tan ^[(9n /35)/(l + »cy )] 
s s 

- (tt/ 2) + C{[2y2g/a + - D (7.53) 

By using Eqs. (7.51) to (7.53), the shock relations [Eqs. (3 ,4) <-(3 .9) ] 

can be expressed, in terms of e3q>anded variables as 

^o- • V - 1/p - (7.54) 

s s s s 

+ Yj^g) -11 (7.55) 


» 3 - ■ - '^^23/'^ - ^ 

Pg- - Pg+ + (1 - l/Pg.) + (1 - I/P3-) 

• ti - 27^3/ a + Yis)!"} 

'* 3 - ■ *‘ 3 + * 

In Eqs. (7,54) through (7,58), only and 
first terms of their expansion. Thus, a series 


(7.56) 

(7.57) 

(7.58) 

u involve y„ in the 
s •*2 

s 

e:q)ansion for the flow 
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varieibles Is assumed al)Out the aucis of symmetry with respect to non* 
dimensional distance ^ near the stagnation streamline as 


p(5,n) = p^ Cn) + P2(h)^^ + — 

C7,59a) 

u(5,n) - u^ (n’)? + — 

(7.59b) 

v(5,n) = v(n) + — 

(7.59c) 

p (^,n) = (n) + — 

(7.59d) 

T(c,n) = T^ (n) + — 

(7.59e) 

u (5,n) = (n) + — 

(7.59f) 

K(?,n) = (n) + — 

(7.59g) 

c (C,n) = c , (n) + — 

P pi 

C7.59h) 

c^(C,n) = (n) + — 

(7.59i) 


Since is a function of downstre 2 ua flow, it cannot be determined by 
2s 

the stagnation solutions. Thus, a value of y^^ * 0 is assumed initially. 
This assumption is removed by iterating on the solution by using the 
previous shock standoff distances to define y 2 g* 

The new form of x-momentum and energy equations in the C,n can be 
written as 


a^w . 3W ^ ^ 

+ a. ^ + a W + a « 0 

37,2 1 a'l 2 3 


(7.60) 


For x-momentTom, W = u and coefficients a^, a^, and a^ 2 ire exactly the 
same as given in [15] , For the enthalpy equation, W = H and again a 


and a^ are the same as defined in [15] but a^ is given by 


+ 2V(1 + + ^yis^r^s’^s^l'^^^^s^”s^ ^^^1'^^^^ 


(7.61) 

Other stagnation streamline equations are the same as given in [15] . It 


should be noted that at the body surface p^ = 1 and p^ = 0. 
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Pig. 6 Finite difference representation of flow field. 



As mentioned earlier, the governing second'-order partial different* 
ial equations are solved by employing an implicit finite-difference 
method. For this purpose, the shock layer is considered as a network 
of nodal points with a veuriaUble grid space in the n'*direction . The 
scheme is shown in Fig. 6 vriiere m is a station measured along the body 
surface and n denotes the station normal to the body surface. The deriv- 
atives are converted to finite-difference form by using Taylor's series 
expeuisions. Thus, unequal space central difference equations in the 
n -direction at point m, n can be written as 


9W 


An, 


n-l 


An. 


w 


w 


3n n An (An , + An ) m,n+l An , (An , + An ) m,n-l 
n n-l n n-l n-l n 


An - An , 

+ n n-l ^ 

An An , m,n 
n n-l 


(7.62a) 


i!«, 2 „ . 2 ^ 

2 n An_ (An_ + An_ , ) m,n+l An_An_ , m,n 


9n 


n n n-l 


n n-l 


W 


An , (An„ + An , ) m,n-l 
n— 1 n n— 1 


(7.62b) 


W — W 

9W . _ m,n m-l,n 

9C'm ■ AC 


(7.62c) 


A typical difference equation is obtained by substituting the above 

equations in Eqs. (7.47) or (7.60) as 

W = - (D /B ) - (A /B )W ^ - (C /B )W - (7.63) 

m,n n n n n m,n+l n n m,n-l 

where 


A = (2 + a An„ ,)/[An„(An + An„ .)] 

n 1 n— X n n n— x 
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B - - 12 - a (An„-&n^ )]/(An An ) - a -- , 

n in n— 1 n n— l 2 4 m— 1 


c « (2 - a.An„)/[An . (An„ + An„ ,)] 

n In n~l n n^l 


3 4 m~*lfn tt"*! 


Now, If It Is assigned that 


W - E W ., + P 
m,n n m,n+l n 


(7.64) 


or 


if *• E W + F 

m,n-l n-1 m,n n-X 


(7.65) 

then by substituting Eq. (7.65) into Eq. (7.63), there is obtained 


W - E-A /(B + C E J)W . 
m,n n n n n>-l m,n+l 

+ (-D - C F , )/(B + C E , ) 
n n n-1 n n n-1 

By comparing Eqs. (7.64) euid (7.66), one finds 


E - -A / (B C E _ ) 
n n n — n n-1 


P 


(-D -CP , )/(B + C E , ) 

n n n-1 n n n-1 


(7.66) 

(7.67) 

(7.68) 


Now, since E. and F. are known from the boundary conditions, E and 
11 n 

F cam be calculated from Eqs. (7.67) and (7.68). The quantities W 
n in ,n 

at point m, n cam now be calculated from Eq. (7.64). 

The overall solution procedure starts with evaluation of the flow 
properties Immediately behind the shock by using the Hankine-Hogoniot 
relations. With known shock and body surface conditions, each of the 
second-order paurtial differential equations are integrated numerically 
by using the tridiagonal formalism of Bq. (7.47) and following the pro- 
cedure described by Eqs. (7.63) to (7.68). As mentioned before, the 
solutions are obtained first for the stagnation streamline. With this 
solution providing the initial conditions, the solution is marched down- 
stream to the desired body location. The first solution pass provides 
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only an approximate flovr-fleld solution. This is because in the first 

solution pass the thin shock<-layer form of the normal momentum equation 

is used, the stagnation streamline solution is assumed to be independent 

of downstream influence, the term dy /d^ is equated to zero at each body 

s 

station, auid the shock angle a is assumed to be the same as the body 
angle 8. All these assumptions eure removed by making additional solu- 
tion passes . 

In the first solution pass, the viscous shock-layer equations are 
solved at any location m after obtaining the shock conditions (to estab- 
lish the outer boundary conditions) from the precursor region solutions. 
The converged solutions at station m-1 are used as the initial guess for 
the solutions at station m. The solution is then iterated locally until 
convergence is achieved. 

For the stagnation streamline, guess values for dependent varieJsles 
are used to start the solution. In the first local iteration, both 
Oy^/8C) and (3W/d^) eire assumed to be zero. The energy equation then 
is integrated numerically to obtain a new temperature. By using this 
temperature, new values of thermodynamic and transport properties are 
calculated. Next, the x-momentum equation is integrated to find the u 
component of velocity. The continuity equation is used to obtain both 
the shock standoff dist 2 mce and the v component of velocity. The 
pressure P is determined by integrating the normal momentum equation. 

The equation of state is used to determine the density. For example, 
the integration of the stagnation streamline continuity equation from 
0 to Ti results in 


[(i + Yi^n) 


D V P. 1 V 

^Is IsU^ 1 




( 7 . 69 ) 
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where 


A » / a + 

0 


C7.70) 


This equation gives the v-velocity component along the stagnation 
streaualine. However, integration of the continuity equation from 
T) » 0 to n * 1 results in 

where 


.1 , 


.1 


B - i p^udn, c - ^ Pj^u^Tidn 

The shock standoff distance can be obtained from the solution of 
Bq, (7.71) as 


Is 


-(2Vi^f2Bu^^)-K(2Vi^^-2BUj^s>^’4Cv^^H-2Cu^^lv^^] 

2 (v +2Cu ) 

Is Is 


1/2 


(7.72) 


Similarly, other quantities at the stagnation streamline are obtained. 

With known stagnation streamline solution and body surface and 
shock conditions, the aibove procedure is used to find solutions for any 
body location m. The downstream shock standoff distance and the v-velocity 
component au:e obtained by integrating the continuity equation in the 
n-direction from 0 to 1, and o to n respectively. Integration of the 
continuity equation from n - 0 to n * 1 results in 

cos 0 p u f pundn + y rp u f pudn 
95 ■'s s s o -"s s s o 


■ (r + y cos 0) [yip u - (1 + y k)p v 1 (7,73) 

s sss sss 

By defining, for station m 

1 — 1 — 

C, “ cos 0 p u / puTidn, C_ = rp u f pudn 
1 ^sso '2 sso 

and denoting the same relations by and for station m-1, Eq, (7,73) 
can be expressed in terms of a difference equation as 
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I(Ciy2 + - (Cjyl + V3)„_iVA5 

■ rp u y ' + cos 8puy*y -rpv 

s s‘'sm 3 s sm sm s s 

- rp V Ky - cos 0 p v y - cos 9 p v icy^ 

sssm sssm sssm 

This C 2 U 1 be expressed in a quadratic form as 

(AA)y^ + (BB)y + (CC) » 0 
sni sni 


where 


AA ■ C, + cos 0Kp V AC 
1 s s 


(7.74) 


(7,75) 


BB s» C- + rp V kAC - cos 0 p u y'AC + cos 0 p v AC 
2ss sss ss 

CC . -tC3(y^)2_^ + 

The shock standoff distance at station m is obtained from Eg. (7.75) as 

y * { - (BB) + C(BB)2 - 4(AA) (CC)]^'^^}/2(AA) (7.76) 

sm 

The V'-velocity component ceui be obtained in a similar manner, Integra* 
tion of the continuity equation from 0 to y gives 

3 7) 

TT 0)p u pudni 

0 sm sm s s 


+ (r + y n cos 9) [(1 + ny k) (p v pv) 
sm sm s s 

“ 0 (7.77) 

sm s s 

As before, this can be expressed in terms of a difference equation as 

{[(KK) - (KK) ,]/AC> + (II) V + (JJ) « 0 (7.78) 

m m*i m m 

where 

(II)^ = (r + y n cos 9) (1 + y^„nic)p v p 
m sm sm s s 

(JJ)^ = -(r + y n cos 0)y’ np„u pu 
m sm sm s s 

(KK)^ ** + y^^n cos 0)p u pudn 

m 0 sm sm s s 

Thus, the v-velocity component at each point on the station m cam be 
obtained from Eq. (7.78). Other quantities at station m are obtained 
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by « similar mannsr. As mentioned before, the first pass Is only an 
approximate solution because of several inhcnrent assumptions. These 
assumptions are resKsved by iteration in the next pass. For the sub- 
sequent solution passes, the shock angle and y are given by 

o - 0 + tan"‘^Iy' /(I + »cy ).J C7.79) 

sm sm 

y ,2 “ C7.80) 

The flow diagrams for computation procedure are shown in Fig. 4 and 
Pigs. 7 to 13. 
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Fig. 7 Flow chart for subroutine SHOCK for shock-layer 
solution. 
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Fig. 9 Flow chaxt for subroutine ENERGY for shock-layer 
solution . 
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Fig- 10 Flow chart for subroutine MOMENTM for shock-layer 
solution. 
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Fig. 11 Flow chart for subroutine RADIATION for shock- layer 
solution. 
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Fig. 12 Flow chart for sxibroutine CHEMIST for shock-layer 
solution . 
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Fig. 13 Definition of integrals used in subroutine RADIATION. 
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8. RESULTS AND DISCDSSXON 


For this study, the entry body considered is a 45^ hyperboloid 
blunt body which enters the Jovian atmosphere at a zero degree angle of 
attack. The body surface is assumed to be gray having a surface emit* 
tance of Q.8. Unless specified otherwise, the surface temperature is 
taken to be imiform at 4,564K. For the case of chemical equilibrium in 
the shock layer, all results were obtadned by considering a body nose 
radius of R* « 23 cm. For chemical nonequilibrium conditions, however, 
three different nose radii (12, 23, and 45 cm) were considered. The 
ncninal composition of the free stream atmosphere was considered to be 
85 percent hydrogen and 15 percent helium for most calculations. However, 
compaurative results were also obtained for the 89' percent hydrogen and 
11 percent helium nominal atmosphere. 

First, results au:e presented for Veuriation in flow properties only 
in the precursor region. These results were obtained with known values 
of radiative heat flux at the shock front. Next, chemical equilibrium 
shock layer results, obtained by considering slip conditions, are pre- 
sented. With these results providing the basis for further investigation, 
complete precursor region-shock layer coupled solutions were obtained for 
both chemical equilibrium cuid chemical nonequilibrium in the shock layer. 
These results acre presented in the last two sections of this chapter. 

8.1 Precursor Region 

In the precursor region, the results were obtained only for the range 
of entry velocities for which free-stream amd shock conditions were avail- 
able (see Table 1). As mentioned before, precursor-region results were 
obtained by employing both the small perturbation method and the thin 
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layer approximation. First small perturbation results are presented, 
and then some key results of this method are compaured with the results 
of the thin layer approximation. 

By employing the small perturbation method, the perturbation quanr- 

titles V- f P- , C - C , and T were calculated numerically amd the 
IZ X H X 

results are illustrated in Figs. 14-^22. In Figs. 14->18, perturbation 
quaintities are shown as a fiinction of distauice fron the shock for differ- 
ent altitudes amd a constant entry velocity of 35 lon/sec. In Figs. 19-23, 
the perturbation quauitities (just ahead of the bow shock) are illustrated 
as a function of the free-stream velocities. Since = “^Iz' separate 
results were not illustrated for the density perturbation. From these 
figures it is evident that the magnitude of perturbation quantities , in 
general, depend on the distamce from the shock, altitude of entry, and 
entry speeds. 

Figures 14-18 show that at a fixed entry velocity, the perturbation 
effects aure greater for lower altitudes and at locations just ahead of 
the shock. This, however, would be expected because the number densities 
of participating species aure greater at lower altitudes and at these 
altitudes most radiative energy from the shock gets absorbed in the 
immediate vicinity of the shock front. At higher altitudes, pertiirba-? 
tion effects are significamt to a laurger distance frcan the shock front. 
This is because, at these altitudes, the number densities of participating 
species au:e small emd radiation effects are felt farther into the free^ 
stream. Specific results presented in Figs. 14-18 indicate that the use 
of the small perturbation theory is justified in determining the velocity, 
density, mass fraction and total enthalpy variations. For example, just 


84 




DISTANCE FROM SHOCK (z/R.) 


Fig. 14 Velocity pertxarbation as a fxanction of distauice 
from the shock at different altitudes euid a 
constant free-stream velocity. 





Fig. 15 Pressiire perturbation as a function of distance 
from the shock at different altitudes and a 
constcuit free-stream velocity. 
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DISTANCE FROM SHOCK (z/F^) 


Fig. 16 Mass fraction of H as a function of distance from 
the shock at different altitudes and a constant 
free-stream velocity. 
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DISTANCE FROM SHOCK {z/R^ 


Fig. 17 Mass fraction of as a function of distance from 
the shock at different altitudes and a constant free 
stream velocity. 





DISTANCE 


Fig. 18 Temperature perturbation s 
the shock at different alt 
stream velocity. 


V. » 35 km /sec 
T « T, + T, 

T« » I45« K 





aQiead of the shock, the value of (V/V ) is 0.9992 for Z « 95 km ^uld Is 

00 

.3 

equed. to 0.99975 for Z - 150 km. Similarly, » 6.8 x 10 for 

-3 ^ 

Z a 95 km and » 2.4 x 10 for Z 150 km (i.e«, 0.68% increzise in 

total enthalpy at 95 km euid 0.24% increase at 150 km) . The static 

pressure and temperature variations, however, cannot be considered small. 

This is because for Z » 95 km, * 2 and = 300 K, and for Z - 150 km, 

» 0.64 auid » 94 K. For these variations, therefore, one could 

q[uestion the validity of the small perturbation theory. 

For different altitudes of entry, perturbation results (just adiead 

of the shock) are illustrated in Figs. 19->23 as a function of entry 

velocities. These results again indicate that the perturbation effects 

are greater for lower altitudes. As would be expected, for any specific 

altitude, the effects are larger for higher entry velocities. This is 

a direct consequence of greater radiative energy tramsfer from the shock 

to the free-stream at high entry speeds. For the most part, variations 

in the velocity, mass fractions, and total enthalpy again are seen to be 

small. For example, for an entry body at an altitude of 95 km, the total 

enthalpy of the gas (H ) entering the shock wave is increased from cdsout 

^1 

0.68 percent at V * 36 km/sec to 1 percent at V = 38 km/sec. For 

Z “ 150 km, however, H increases from 0.24 percent at 35 km/sec to 

^1 

0.66 percent at 42 km/sec. The variations in the static pressure and 
t^perature, in some cases, are seen to be several times greater them the 
ambient values. These Icurge variations, however, occur for conditions 
where dissociation is high and the validity of the entire theory is 
questionable [1,2]. 

By employing the governing equations (£qs. 2.9-2.13) and the spec- 
tral infoannation of Sec. 6.2.1, numerical results were obtained for 
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ALTITUDE (Z 
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FREE-STREAM VELOCITY (V.), km/sec 

Fig. 19 Pressure perturbation (just ahead of the shock) as 
free-stream velocity for constant altitudes. 





FREE-STREAM VELOCITY km/sec 


20 Mass fraction, of H (just ahead of the shock) as a 
function of free-stream velocity for constant 
altitudes . 
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FREE-STREAM VELOCITY (V;,), km/sec 


Fig. 21 Mass fraction of (just ahead of the shock) as a 

function of free-streeun velocity for constant 
altitudes . 
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FREE-STREAM VELOCITY (VJ, km/sec 

Fig. 22 Temperature perturbation (just ahead of the shock) as a functl 
free-stream velocity for constcint altitudes. 










ALTITUDE (Z) • 116 km 

'THIN LAYER 

perturbation 

-V» » 39,09 km sec 

/- 35.00 
' ^30.00 


distance from shock, cm 


Fig. 26 Con^)arlson of results for pressure variation in the 
precursor zone. 
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velocity, presstire and temperature variations for different values of 
n at s * 0, Specific results for an altitude of Z ■* 116 km are com- 
pared in Figs. 24-26 with corresponding results of the small perturba-r 
tion theory. For the reinge of parameters considered, the results 
obtained by the two proced\ires are seen to be in excellent agreement. 

It is obvious from these results that either approach could be utilized 
in the investigation of the precursor region flow field. It was noted 
in Sec, 7.1 that for the Jupiter's entry conditions, the general govern- 
ing equations of the small perturbation theory reduced to the case of 
simple plane source. As such, use of this method to Jupiter's entry 
case is restricted to one-dimensional analyses . The advantage of thin 
layer approximation procedure is that it is physically more convincing 
euid it ceui be extended easily to three-dimensional and cixisymmetric 
cases . 

8.2 Effects of Shock and Body Slip Conditions 
By invoking the boyd auid shock slip conditions , results for Vciria- 
tion in the shock layer flow properties were calculated for higher alti- 
tude entry conditions. Seme important results of this investigation are 
presented in this section. Results are presented first for the velocity 
emd temperature jumps at the body surface. Following this, results are 
presented for the properties immediately behind the shock. Next, the 
effects of radiation on convective heating at higher altitudes cure dis- 
cussed. Finally, to assess the influence of slip conditions, results 
zure presented for the convective and radiative heating. It should be 
emphasized here that the term slip conditions (or slip boundziry condi- 
tions] , as used in this study, implies both the body and the shock slip 
conditions . 
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The variation in the surface slip velocity is illustrated in Fig. 27 
as a function of the entry altitudes. Since u-velocity is almost zero 
at the stagnation streamline, the results presented in Fig. 27 have been 
obtained for loQation (or station) 3 of Fig. 6. Figure 27 cleaurly illus* 
trates that the condition of no slip is not satisfied at higher altitudes. 
Since u is normalized by the shock value (i.e., u » u/u^) , the magnitude 
of velocity slip can be expressed as a percent of u^. It is evident from 
Fig. 27 that about 8 percent velocity slip occurs at Z - 261 amd only 
0.1 percent at Z * 143 km. 

The temperature jump at the body surface is shown in Fig. 28 for 
different entry altitudes. The results presented in this figure are for 
the case with no radiation and, in obtaining these results, the body 
surface temperature was taken to be 4,000 K. A temperature jump of about 
18 percent (i.e., VT = 680 K) is noted at the stagnation point for entry 
conditions at z =261 km. At lower altitudes, however, the temperature 
jtjmp is seen to be relatively small. For example, at Z = 116 km, the 
temperature jump is only 3 K. 

Figures 29-31 show the temperature jump, velocity slip, density cuid 
total enthalpy changes just behind the shock. It is evident from Fig. 29 
that v^en the altitude is lower them 225 km, the shock slip conditions 
are not important. However, a significamt temperature difference is 
noted at Z = 261 km. The results presented in Fig. 30 illustrate that 
both the u and v velocity components are influenced by the slip condi- 
tions. Since both the temperature and velocity components decrease just 
behind the shock, the slip conditions result in an increase in density 
and a decrease in total enthalpy. This is cleaurly evident from the 
result of Fig. 31. 
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Pig. 27 Velocity slip along the body surface for different 
entry ailtitudes . 
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Fig. 28 Temperature jump along the body s\orface for 
different entry altitudes. 
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Figure 32 shows how the convective heat flux is effected by radia- 
tion. The results indicate that at low altitudes, the convective heat 
flux decreases with increasing altitude and at high altitudes it increases 
with the altitude. This is because at different ranges of altitude, the 
temperature distribution is rearranged by the radiation effects. It is 
noted that a maximum of 50 percent change in convective heat treuisfer 
occurs at Z » 261 km auid a 25 percent change at Z = 225 km. 

Figure 33 shows how the radiative heat flux is affected by slip 
boundary conditions. It is seen that the effect is very small at alti- 
tudes lower them 225 km. It is foxmd that there is approximately 50 
percent reduction in radiative heat flux due to the shock temperature 
jump at 261 km. 

The effects of slip boundary conditions on convective heat flux 
towards the body (along the body surface) are illustrated in Figs. 34 
emd 35, for the cases with and without the radiation interaction. The 
results indicate that the slip conditions start to effect the convective 
heat flux at Z = 225 km (e = 0.09064) by approximately 8 percent (at 
stagnation point) and this increases to 27 percent at Z = 261 km 
(e = 0.2129). The effect is seen to increase with the distance away 
from the stagnation point, and inclusion of radiation is seen to suppress 
this influence. When the altitude is less than 225 km, the effect of 
slip boundary conditions is relatively small aind it Ccin be neglected. 

8.3 Influence of Precursor Heating 
on Viscous Equilibrium Flow 

By considering the conditions of chemical equilibrium in the shock 
layer, governing equations of both the precursor and shock layer regions 
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were solved for physically realistic Jovian entry conditions. Results 
of complete parametric study eure presented in this section. First, the 
results are presented for quantities just behind the shock wave, 2 uid 
then a few results of flow varietbles within the shock layer eu:e pre- 
sented. Next, results are presented for the entire shock-precursor 
region. Finally, a few results are presented to demonstrate the influ- 
ence of precursor heating on the magnitude of different heat fluxes in 
the shock layer. 

The radiative flux from the shock layer towards the precursor 
region is found to be highest at the stagnation line shock location. 
Results of the radiative flux from the shock front are shown in Fig, 36 
for different altitudes of entry. As would be expected, prec\irsor 
heating results in a higher radiative flux at the shock front. It is 
seen that the radiative flux reaches a maximum value for an altitude 
of adsout 116 km, and the largest precursor effect (PE) of about 8 percent 
is found to be for this altitude. This is a direct consequence of the 
free strecun emd entry conditions at this altitude. For other entry 
conditions (altitudes ) , precursor effects are seen to be relatively 
lower. 

Figxire 37 shows the shock stcindoff Vciriation with distance along the 
body surface for different entry altitudes. The shock standoff distance, 
in general , is seen to decrease with increasing altitudes . This is 
because higher entry velocities aore associated with higher altitudes. 

The precursor heating results in a slight increase in the shock standoff 
distance (a maximum of about 2 percent for Z = 116 km) because the 
density of the shock layer is slightly reduced. 


Ill 






5-x7 r* 


Fig. 37 Shock stauidoff distiance veuriation with distance along 
body surface. 
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The conditions jusrt behind the shock aixe illustrated in Figs. 38^1 
as a function of distsmce along the -body for different entry altitudes. 
For z » 116 km. Fig. 38 shows that precursor heating increases the 
enthalpy by a Tnayimum of about 2 percent at the stagnation line. The 
chemge in shock temperature is shown in Fig. 39 for different altitudes. 
As would be expected, precursor heating results in a relatively higher 
temperature. The effect of precursor heating on the pressure just behind 
the shock was found to be small and, therefore, it could not be shown in 
a figure conveniently. Since the pressiire essentially remains unchanged, 
precursor heating results in a decrease in the density (see Fig. 40) 
mainly because of an increase in the temperature. It was found that 
precursor heating had no significeuit influence on the u-oompohent of 
velocity, but the v-component is slightly increased (see Fig. 41) as a 
result of decrease in the shock density. 

Variations in pressure, density, velocity, and. chemical species 
across the shock layer are shown in Figs. 42-^4 for cui altitude of 
Z « 116 km. Results presented in these figures are normalized by their 
shock values and they show that precursor effects are felt throughout 
the shock layer. Results presented in Figs. 42 and 43 for two body 
locations (C = 0 and 1) indicate the relative change in pressure, density, 
and velocities as compared to their shock values. For 5=0, Fig, 44 
shows that precursor heating slightly decreases the concentration (mole 
fraction) of atomic hydrogen and increases the concentration of ions and 
electrons throughout the shock layer. 

Variations of temperature, pressure, density, and velocity along 
the stagnation streamline in the entire shock layer-precursor zone are 
illustrated in Figs. 45-48 for different altitudes. Since higher entry 
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Pig. 38 Enthalpy variation just behind the shock with distance 
along the body surface. 
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Fig. 41 Variation of v-velocity component just behind the 
sh<~^k with distance along the body surface. 
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Fig. 42 V 2 u:iation of pressure and density in the shock layer 
for two body locations (5*0 auid 1) . 
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Fig. 43 Variations of velocity co■^ponent8 in the shock layer for 
locations (C ** 0 and 1) . 
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Fig* 44 Species concentration in the shock layer for 
S =* 0. 


121 


r. K 

6 500{ 


PRECURSOR EFFECT 

NO PRECURSOR EFFECT 


2 = 103 km 





•SHOCK 


PRECURSOR ZONE 


z = 143 km 



LO 3L0 5.0 T.OaO 

n = y/yj = y*fy* 


Fig, 45 Temperature variation in the shock/precursor 
region along the stangation streamline. 
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Pig. 47 Density variation in the shock/precursor region 
along the stagnation streamline. 



Pig. 48 Variation of v- velocity con^jonent in the shock/ 
precursor region along the stagnation streamline. 
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velocities are associated with higher altitudes , precursor effects , in 
general, are found to be larger for higher altitudes. The results for 
the precursor region show a dramatic increase in the pressure and 
temperature but only a slight change in the density and velocity. The 
changes are largest near the shock front because a major portion of 
radiation from the shock layer gets absorbed in the immediate vicinity 
of the shock front. Figures 45 and 46 show that, in spite of a large 
increase in the temperature and pressure in the precursor region, pre*?' 
cursor heating does not change the temperature and pressure distribution 
in the shock layer drama ticad.ly. The change in temperature, however, is 
significant and Cas would be expected) the maximum chcuige occurs just 
behind the shock. There is a slight change in the pressure near the body 
but virtually no change closer to the shock. Figure 47 shows that the 
change in density in the shock layer is higher for higher altitudes and 
towards the shock. As discussed before, precursor heating results in a 
slight decrease in the shock layer density. Virtually no change in the 
u-component of shock layer velocity was found, but, as shown in Fig. 48, 
the v-component is slightly increased. 

The effects of precursor heating on different heat fluxes in the 
shock layer are illustrated in Figs. 49-51. These results clearly 
demonstrate that precursor heating has a significant influence on 
increasing the heat transfer to the entry body. This increase essentially 
is a direct consequence of higher shock layer temperatures resulting from 
the upstream absorption of radiation. Figure 49 shows the variation of 
radiative and convective heat flux with distance along the body surface 
for Z = 116 km. It is noted that the precursor heating results in a 
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Fig. 49 Vsuriation of radiative ^md convective heat flxix with 
distance along the body surface for Z = 116 km. 
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Fig. 50 Variation of radiative heat flux in the shock-layer 
for two body locations (C * 0 and 0.5), 2 « 116 km. 
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Fig. 51 Radiative auid convective heat flux to the body (at 
5*0, n » 0) for different entry altitudes. 



7.5 percent Increase in the radiative flux and adsbut 3 percent increase 
in the convective flux to the body at the stagnation point. The increase 
in heat transfer at other body locations axe relatively lower. A 
similar conclusion can be dravm from the results presented in Fig. 50 
for the radiative flux towaurds the shock auid the body for two body 
locations (S “ 0 and 1) at Z » 116 km. Results of radiative amd con-<* 
vective heat flux at the body (for ? « 0) aure illustrated in Fig, 51 for 
different altitudes of entry. The radiative flux results aire seen to 
follow the trend exhibited in Fig. 36 for radiation at the shock front. 
The convective heat flux, however, is seen to increase slowly with the 
altitude up to Z = 131 km and thereafter decrease with increasing alti- 
tudes. The precursor effect is found to increase the radiative heating 
by a maximum of aibout 7.5 percent at Z ^ 116 km auid the convective 
heating by 4.5 percent at Z =» 131 km. 

8.4 Influence of Precursor Heating 
on Viscous Nonequilibri\am Flow 

The influence of precursor heating on the flow phenomena euround a 
Jovian entry body was investigated under the conditions of chemical non-^ 
equilibrium in the shock layer. As mentioned before, the entry body 
considered for this study is a 45® hyperboloid blunt body. The body 
enters the Jupitor's atmosphere at zero angle of attack. The two nominal 
atmospheres considered for Jovian entry consist of 85 and 89 percent hydro- 
gen (by mole fraction) respectively. Also, to investigate the influence 
of change in the body nose radius on the thickness of the nonequilibrixjm 
layer, three different nose radii (12, 23, and 45 cm) were considered. 

To illustrate the importeuit features of the nonequilibrium euialysis, most 


130 



results were obtained for entry conditions which closely correspond to 
the peak heating conditions (i*e., for conditions at Z “ 116 km). How- 
ever, a few illustrative results have also been obtained for other 
entry conditions. Equilibrium and nonequilibrium results are presented 
first for variation of different properties in the shock layer. Results 
are then presented to illustrate the influence of precursor heating. 
Finally, results are presented for variation of different heat fluxes 
in the shock layer under the influence of both the nonequilibrium con- 
ditions and the precursor heating. 

Two assiomptions can be made 2 ibout the moleculeu: hydrogen entering 
the shock layer immediately behind the shock. One criteria is to assimie 
that chemical reactions are "completely frozen" and initial composition 
of hydrogen just bePiind the shock corresponds to the free stream value. 
The second criteria is to consider that all hydrogen molecules have been 
dissociated immediately behind the shock. This is referred to as the 
"half o frozen” condition. Nonequilibrium results obtained for these two 
cases Cfor entry conditions at Z * 116 km emd for 85 percent hydrogen 
nominal atmosphere} are illustrated in Figs. 52 to 54 as a function of 
the normal coordinate at the stagnation point. Figure 52 shows the mole 
concentration of different species across the shock layer. Xt is evi- 
dent from this figure that molecular hydrogen is completely dissociated 
within 2 Lbout 4 percent of the total shock stauidoff distance from the 
shock wave. This is referred to as the dissociation zone (or the 
dissociated region) , The Veiriation in nondimens ional v-velocity compon- 
ent euid density is illustrated in Fig. 53. Since molecular weights 
change rapidly in the dissociated region, there is ^ increase in 
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Fig. 53 Equlllbriisa and noneguilibrlum v-velocity con^onan't 
amd density variation in the shock layer for ^ * 0^ 







velocity and a decrease in density near the outer edge of the dissocla->> 
tlon zone . The temperature distribution Is shown In Fig . 54 . It Is 
noted that the temperature just behind the shock wave reaches a value 
of approximately 45,000 K in the completely frozen condition. After a 
short Interval, however, eU.1 hydrogen molecules acre dissociated and 
temperature drops to about 25,000 K. Next, Ionization occurs and, as a 
result of this, temperature continues to decrease until it reaches the 
equilibrium value. From the results presented In Figs. 52^54, it Is con<^ 
eluded that the half' frozen and completely frozen assumptions are quite 
close except in the dissociated region neaur the shock wave, and that 
the half-frozen flow computation is a reasonably good assun^tlon for 
conditions of chemical nonequilibrium at altitudes near the peak heating 
region. Thus, all other results presented in this section have been 
obtained by considering only the half- frozen condition behind the shock. 

As discussed in the previous section, the shock stauidoff distauice 
Cfor a given body nose radius) varies with the altitude of entry and 
entry velocity. It should be pointed out here that, in general, the 
shock standoff distance increases with increasing the body nose radius . 
For entry conditions at Z » 116 km, equilibrium and nonequilibrium 
results for the shock standoff distance are illustrated in Fig. 55 as 
a function of the coordinate along the body surface. It is noted that 
the shock standoff distances for equilibrium and with radiation axe con^ 
sideredsly lower than for nonequilibrium and with no radiation. This, 
however, would be expected because shock-layer densities are greater for 
radiation and equilibrium conditions them for no radiation and nonequlll- 
brium conditions . 
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Fig. 55 Shock standoff variation with distance along the body 

surface. 
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Variations in chemical species across the shock layer are shown 
in Figs. 56-^59 for different conditions. For entry conditions at 
Z ■ 116 km, results presented in Figs. 56 and 57 show that the non<r> 
equilibrium layer is about 25 percent of the total shock- layer thickness 
for no radiation case and about 50 percent for the case with radiation. 
This is because inclusion of radiation results in a different tempera- 
ture distribution in the shock layer. This point will be discussed 
further %diile presenting results for the ten^erature variation. Near 
the wall, the mass fractions of atomic hydrogen and electrons eu:e higher 
for nonequilibrium conditions with radiation, lliis is because cold gases 
near the wall £dbsorb relatively more radiative heat flux in nonequili- 
brium case. For no radiation case, a comparison of results presented 
in Figs. 56 amd 58 reveal that the noneqiiilibrium layer increases from 
25 percent at Z « 116 km to aibout 40 percent at Z > 143 km. This is 
because density is lower at higher altitudes and, therefore, it will 
take a relatively longer time to reach equilibrium condition. For 
Z ■ 116 km entry conditions. Fig, 59 shows the species concentrations 
for three different body nose radii (12 , 23 , and 45 cm) . These results 
Indicate that the thickness Cor range) of the nonequilibrium layer de- 
creases with increasing nose radius. In particular, it is seen that the 
thiclcnesa is about 40 percent for 12 cm but it is only 10 percent 

for * 43 cm. This is because the shock stauidoff distauice is propor- 
tional to the body nose radius emd the releucation time for chemical 
reactions is zdx^ut the same for all cases. 

Temperature distributions across the shock layer sure illustrated in 
Figs. 60-62 for different conditions. For the case with no radiation, 
the heavy particle auid electron temperature variations across the shock 
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Fig. 56 Species concentration in the shock layer for 
C “ 0 (with no radiation) . 
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Fig. 57 Species concentration in the shock layer for K 
(with radiation) . 









tlon} in Pig. 60 for different entry conditions. The results show that 
in the 2 Lhsence of radiation, the noneguilibrium temperature is higher 
than the equilibrium temperature throughout the shock layer for each 
entry condition. It is also noted that the electron temperature, which 
is lower theui the heavy pcurticle temperature during early stages of 
ionization, asymptoticeLLly approaches the heavy particle temperature 
during the later stages of ionization. As discussed in the previous 
section, the temperature distribution in the shock layer is relatively 
higher for higher altitudes because of higher entry velocities. For 
entry conditions at Z - 116 km, the electron temperature distributions 
(without and with radiation) are shown in Figs, 61 and 62 for three 
different body nose radii. As noted earlier, the thickness of the non--> 
equilibrium layer decreases with increasing nose radius. Also, for a 
given nose radius, inclusion of radiation increases the thickness of the 
nonequilibrium layer. This is because the loss of radiation from the 
shock layer results in an entirely different temperature distribution 
(see Fig. 62) amd leaves relatively less energy for dissociation amd 
ionization of the gas. 

For entry conditions at Z => 116 km. Fig. 63 shows the mass fraction 
of atomic hydrogen auid hydrogen ion along the stagnation streamline in the 
precursor region. While equilibrixmi results indicate that only 5 percent 
hydrogen is dissociated amd 0.018 percent is ionized, the nonequilibrium 
results show that 15 percent hydrogen is dissociated and 0.8 percent 
ionized. It should be pointed out that the composition of the precursor 
gas will be different for different entry conditions. It should be 
emphasized here again that in investigating the precursor region flow 
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Fig. 61 Temperature variation for different body nose 
radius in the nonequilibrium region at ^ « 0 
(with no radiation) . 
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Fig. 62 Temperature vairiation for different body nose radius 

in the nonequilibrium region at C * 0 (with radiation) 
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properties and their influence on the shock layer flour phenonena^ the 
entire precursor-i-shock layer solutions are obtained by iterative protr 
cedures . 

For the case with radiation euid for entry conditions at Z « 116 taiy 
the heavy particle and electron temperature variations across the shock 
layer are illustrated in Fig. 64 along with the equilibrium temperature 
distribution. In ccmpsurisic^ with results of Fig* 60, it is seeh that 
in the present case, the nonequilibriTsa t^perature is lower than the 
equilibrium temperature in certain portions of the shock. This is a 
direct consequence of the radiation cooling (i . e . , radiation loss to the 
free stream) of the shock layer. Also, in this case, the nonequilibriun 
temperat\u:e is slightly higher than the equilibrium temperature in the 
vicinity of the wall. This is because cold gases near the wall absorb 
radiation from the high temperature region of the shock layer. As would 
be esqiected, precursor heating results in a slightly higher shock-layer 
temperature distribution. 

Variations of temperature, pressure, and density along the stagna** 
tion streamline in the entire shock layer-precursor zone aure illustrated 
in Figs. 65-68 for different conditions. These results show that pre- 
cursor effects are higher for the nonequilibrium conditions. This, 
however, would be expected since in this case, the radiative heat flux 
towards the precursor region is considerably higher. The shock- layer 
nonequilibrium condition significantly influences the temperature and 
pressure vairiations in the precursor zone, but its effects on density 
cheuiges axe quite small. As noted earlier, in the shock layer, non- 
equilibriimi results approach the corresponding equilibriimi values at 
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Fig. 64 Equilibrium and nonequilibrium temperature variation 
in the shock layer (with radiation? . 
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Fig. 65 Temperature variation in the shock/precursor region 
2 d.ong the stagnation streamline, Z « 116 km. 
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Fig. 66 Temperature variation in the shock/precursor region 
along the stagnation streamline, Z = 143 km. 
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Fig. 67 Pressure variation in the shock/precursor region 
along the stagnation streamline. 
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p = p */p« 



Pig, 68 Density variation in the shock/precursor region along 
the stagnation streamline. 
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about 25 percent of the shock, layer thickness frcm the shock wave. For 
the equilibrium case, the influence of precursor heating on shock- layer 
temperature, pressure, and density veuriations is discussed in. the pre^ 
vious section. 

For a comparision of the shock-layer flow phenomena for the two 
nominal compositions of the Jovian atmosphere, illustrative results were 
obtained for entry conditions at Z « 116 km. Results for the temperature 

s. 

variation immediately behind the shock auid for the radiative heat flux 
across the shock layer eire illustrated in Figs, 69-71. It is evident 
from Fig. 69 that the shock temperature is lower by about 2 percent for 
case of 89 percent hydrogen atmosphere. This is because, in this case, 
relatively more energy is required to dissociate the molecular hydrogen. 
Since the shock temperature is lower in this case, the radiatve heat 
fluxes (q^ as well as q ) are lower for both equilibrium and nonequili- 
brium conditions (see Figs. 70 and 71), 

To investigate the extent of heating on an entry body, the vauria- 
tions in radiative heat flux in the shock layer were calculated for 
different conditions. As discussed earlier, the chemical nonequilibriiam 
effects are more importamt with small body nose radius and for higher 
altitude entry conditions. Results for radiative flux towards the shock 
and body are shown in Fig. 72 for R* = 12 cm and 2 = 116 km. The 
results indicate that, in the nonequilibrium case, the radiative heat 
flux is increased to about 70 percent toward the body and almost 2,5 
times toward the shock (i.e., toward the precursor region). Results for 
radiative heating of the body for R* = 23 cm and Z » 143 km are shown in 
Fig. 73. The results show that the heat flux is about three times higher 
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Fig. 73 Equilibrium and. nonequilibrium radiative heat- flux 
variation in the shock layer, Z = 143 km. 





for the nonequilibrium conditions. This is a direct consequence of the 
higher temperature in the nonequilibritsn layer near the shock. 

To investigate the influence of precursor heating on viscous non-- 

equilibrium shock— layer flow phencmena, specific results were obtained 

for the peak heating entry conditions auid for an entry body with a nose 

radius of R* * 23 cm. These are presented here as final results of the 
n 

present study. 

The radiative heat flux from the shock layer towards the shock front 
and the precursor region is shown in Pig, 74 for both equilibrium auid 
nonequilibrium conditions. The results cleaurly indicate that heat flux 
toward the precursor region is considercibly higher for nonequilibrium 
conditions. This is again a direct consequence of higher nonecjui librium 
temperature in the shock layer. As discussed before, precursor heating 
results in a higher radiative flux at the shock front. The results of 
Fig. 74 indicate that precursor heating results in a 15 percent increase 
in radiative flux in the nonequilibrium case while only 8.5 percent 
increase is noticed for the equilibrium condition. 

The results of equilibrium cind nonequilibrium radiative flux towards 
the body (along the stagnation line) are illustrated in Fig, 75, Although 
it is realistic to calculate the radiative flux based on the electron 
temperature, results (for the case with no precursor effects) have been 
obtained also by using the heavy particle temperature only for compara- 
tive purposes. The nonequilibrium results are seen to be significantly 
higher than the equilibrium results. This, is primarly due to the high 
temperature region near the shock where nonequilibrium temperature over-i- 
shoots occur . 
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Fig. 74 Equilibriiam and nonequilibrium radiative heat flux 
towards the shock for 5 » 0. 
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Flgiire 76 shows the veuriation of radiative convective fliuc with 
distamce along the body surface. The radiative as well as convective 
heat transfer to the body surface is seen to be enhanced by the non^ 
equilibrium conditions. . As discussed above, the increase in radiative 
heating is a direct consequence of higher electronic temperature. For 
the case with no radiation, the convective heat flux toward the body was 
found to be the same for equilibriimi euid nonequilibrium conditions. For 
the case with radiation, however. Fig. 76 shows that the convective heat 
flux for the nonequilibriiam case is about 20 percent higher thain the 
corresponding equilibrium value at the stagnation point. This is because 
the cold gas near the wall absorbs higher radiative flux from the shock 
layer under the nonequilibrium conditions. 7^ discussed before, the 
influence of precursor heating is enhanced due to nonequilibrium condi^ 
tions. Figure 76 shows that precursor heating results in a 10.5 percent 
increase in the radiative flux at the stagnation point in the nonequili- 
brium case while only about 7 percent increase is noted for the equili- 
briiam case , 

For the entry conditions considered in this study, therefore, it 
is logical to conclude that nonequilibrium heating of the body is signi- 
ficantly higher than equilibrium heating. Results similar to this were 
also obtained by Grose and Nealy 176] for Venusian entry conditions. 

For certain Jovian entry conditions, results presented in 18,59] indicate 
that nonequilibrium heating is considerably less than the ecpjilibrium 
heating. This obviously is in contradiction to the present findings. 

It should be pointed out that for the entry conditions considered in 
this study, the temperature just behind the shock is very high and all 
hydrogen molecules aure completed dissociated . Under these conditions , 
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Pig. 76 Variation of radiative and convective heat flux with 
distance along the body surface. 







H line emissions are higher th£ui for the equilibrium conditions. This 
is because, in addition to high temperature, the number density of atomic 
hydrogen is consider cUbly higher thzm the equilibriim value. Thus, find^ 
ings of the present study appear to be completely justified. 


164 



CONCLUDING REMARKS 

The main objective of this study was to investigate the influence 
of precTxrsor heating on the entire shock layer flow phenomena around a 
Jovian entry body under physically realistic conditions. For this pur--' 
pose, the flow in front of the entry body was divided into three regions, 
the shock layer, the precursor zone, auid the free stream. The problem 
was formulated by considering the chemical equilibrium as well as non- 
equilibrium composition of the shock layer gas. 

In the precursor region, flow phenomena was investigated by 
employing the small perturbation theory of classical aerodynamics and 
the thin layer approximations of hypersonic flow. For Jovian entry 
conditions, one -dimensional results obtained by the two methods were 
found to be in good agreement for the remge of parameters considered. 

The results, in general, indicate that for certain combinations of entry 
speeds euid altitudes of entry, the prec\arsor effects cannot be ignored 
while cmalyzing flows around Joviaui entry bodies. The usefulness of the 
thin-layer approximation in analyzing the precursor region flow is 
demonstrated. The main advamtage of this method is that it is physically 
more convincing and its use cam be extended easily to aixisymmetric and 
three-dimensional cases . 

In the shock layer, results of flow variables were obtained along 
the body and the bow shock and across the shock layer for different 
entry conditions. The results show that the slip boundary conditions 
(both at the shock wave amd the body) should be used when the entry 
altitudes aore higher than 225 km. Specific results for the chemical 
equilibrium condition indicate that, in most cases, precursor heating 
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has a maximum influence on flow variables (except the pressure) at the 
stagnation line shock location. It was found that while pressure 
essentially remains unchanged in the shock layer, the precursor heating 
results in cui increase in the enthalpy, temperature, and v-component of 
velocity, and a decrease in the shock layer density. For the entry con- 
ditions considered in this study, results clearly demonstrate that pre- 
cursor heating has a significant influence on increasing the heat 
transfer to the entry body. Chemical nonequilibrium results reveal that 
there exists a noneguilibrium layer of considerable thickness in the 
shock layer region and inclusion of the radiative heat flux term in the 
energy equation increases the thickness of this layer. Under nonequili- 
brium conditions, temperature (heavy paurticle as well as electronic) 
overshoots occur near the shock wave. As a result of this, the radiative 
as well as convective heat transfer to the body surface is increased 
significcintly . The influence of precursor heating is enheuiced due to 
nonequilibrium conditions; a 9.5 percent increase in the stagnation point 
radiative heating has been observed at an altitude of 116 km. 

For further studies, it is suggested that the precursor region flow 
phenomena be investigated without making the tliin layer approximation. 
Since precursor region is relatively thin for most entry conditions, this 
improvement probadsly will not change the findings of the present study. 
However, it would be advisable to consider turbulent flow in the shock 
layer , especially for analyzing the flow away from the stagnation region , 
Also, a more general model for radiative tremsport (instead of the tan- 
gent slab approximation) should be used, and conditions of different 
cuigles of attack for the entry body should be considered. 


166 



REFERENCES 


1. Smithy G, L. , "Radiation- Induced Precursor Plow Field Ahead of a 
Reentering Body," Ph.D. dissertation, Virginia Polytechnic Insti- 
tute, March 1968. 

2. Smith, G. L. , "Radiation Induced Precursor Flow Field Ahead of a 
Reentering Body," AXAA Paper No, 68-667, June 1968. 

3. Liu, J. T. C. , "synoptic: Influence of Upstream Absorption on the 

Stagnation Region shock Layer Radiation," AIAA Journal , Vol. 8, 

No. 10, Oct. 1970, pp. 1730-1737. 

4. Liu, J. T. C. , "Influence of Upstream Absorption on the Stagnation 
Region Shock Layer Radiation," AIAA Journal , Vol. 8, No. 10, 

Oct. 1970, pp. 1730-1737. 

5. Murty, S. S. R. , "Effect of Line Radiation on Precursor Ionization," 
J. Quamtitative Spectroscopy s Radiative Transfer , Vol. 8, No. 1, 
Jan. 1968, pp. 1605-1617. 

6. Nelson, H. F. and Goulard, R. , "Structure of Shock Waves with Non- 
equilibriimi Radiation euid Ionization," The Physics of Fluids , 

Vol. 12, No. 8, Aug. 1969, pp. 1605-1617. 

7. Perri, A. N. and Cl 2 u:ke, J. H. , "Radiative Ionization Patterns in 
Cold Precursor of Axisymmetric Detached Shock , " AIAA Joixcnal , 

Vol. 8, No. 9, Sept. 1970, pp. 1574-1581. 

8. Leibowitz, L. P. and Kuo, T. J, , " lonizational Nonequilibrium 
Heating Diiring Outer Planetary Entries," AIAA Journal , Vol. 14, 

No. 9, Sept. 1976 , pp. 1324-1329. 

9. Livingston, F. R. and Poon, T. Y., "Relaxation Distauice and 
Equilibrium Electron Density Measxarement in Hydrogen-Helium 
Plasmas," AIAA Journal , Vol. 14, No. 9, Sept. 1976, pp. 1335- 
1337. 

10. Garrett, L. B. , "An Implicit Finite Difference Solution to the 
Viscous Radiating Shock Ijayer with Strong Blowing," NASA TMX-67584, 
1971. 

11. l^derson, J. D. , Jr., "An Engineering Survey of Radiating Shock 
Layer," AIAA Journal , Vol. 7, No. 9, Sept. 1969, pp. 1665-1675. 

12. Sutton, K. , "Characteristic of Coupled Nongray Gas Flows with 
Ablation Product Effects About Blxmt Bodies During Planetary 
Entries," Ph.D. Dissertation, North Carolina State University, 

1973. 


167 


13. Kximar, A, ; Tiwari, S. N. ; and Graves, R. A.: Effects of Small 

Angle of Affack on Radiating Viscous Shock Layer Solutions for 
Jovian Entry, AIAA Paper, No. 78-909, May 1978. Also, 

Technical Report NSG<-1464, May 1978, School of Engineering, 

Old Dominion University, Norfolk, Virginia. 

14. Davis, R. T. : Numerical Solution of the Hypersonic Viscous 

Shock-Layer Equations, AIAA Journal , Vol. 8, No. 5, May 1970, 
pp. 843-851. 

15. Moss, J. N. : Reacting Viscous -Shock-Layer Solutions with Multi- 

component Diffusion and Mass Injection. NASA TR R-411, June 1974. 

16. Moss, James N. : Stagnation and Downstrecua Viscous-Shock- Layer 

Solutions with Radiation euid Coupled Ablation Injection. AIAA 
Paper 74-73, AIAA 12th Aerospace Sciences Meeting, Washington, 
D.C., Jem. 1974. 

17. Moss, J. N. : Radiative Viscous-Shock- Layer Solutions with 

Coupled Ablation Injection, AIAA Journal , Vol. 14, No. 9, 

Sept. 1976, pp. 1311-1316. 

18. Lasher, L. E.; and Wilson, K. H. ; Effect of Shock Precursor 
Absorption on Superorbital Entry Heating. AIAA Journal , Vol. 

6, No. 12, Dec. 1968, pp. 2419-2420. 

19. Lasher, L. E.; auid Wilson, K. H. : Effect of Shock Precursor 

Heating on Radiative Flux to Blunt Bodies. NASA CR-1265, 

Feb. 1969. 

20. Hudson, R. D.; Critical Review of Ultraviolet Photoeibsorption 
Cross Section for Molecules of Astrophysical and Aeronomic 
Interest, Review of Geophysics and Space Physics , Vol. 9, No. 

2, May 1971, pp. 305-406. 

21. Lee, L. C. ; Carson, R. W. ; and Judge, D. L.; The /^sorption 

Cross Sections of H„ and D„ from 180 to 780 A, Journal of 

2 2 ' 

Qucuititative Spectroscopy & Radiative Transfer , Vol. 16, No. 

10, Oct. 1976, pp. 873-877. 

22. Nicolet, W. E. ; User's Manual for the Generalized Radiation 
Treuisfer Code (RAD/EQUIL) , Aerotherm Report No. UM-69-9, 

Aerotherm Corp. , Mountain View, California, 1969. Also: User's 

Manual for RAD/EQUIL/1973 , A General Purpose Radiation Trans- 
port Program, NASA CR-132470, 1973, 

23. Wilson, K. H. : RATRAP - A Radiation Transport Code. Report No. 

6-77-67-12, Lockheed Missiles & Space Co., Sunnyvale California, 
March 14, 1967, 

24. Thomas, M.; The Spectral Linear Absorption Coefficient of Gases - 
Computer Program SPECS (H 189), Douglas Report DAC-59135, 
McDonnell -Douglas Astronautics Co., Western Division, Santa 
Monica, California, 1967. 


168 


25. 


Falemga, R. A.*, and Olstad, W. B.: 7\n Approximate Inviscid Radia- 

tion Flow-Field Analysis for Sphere-Cone Venusieui Entry Vehicles, 

AIAA Paper No. 74-758, July 1974. 

26. Zol^, E. V.; Sutton, K. ; Olstad, W. B. ; euid Moss, J. N.* A 
Approximate Inviscid Radiation Flow Field Analysis of Outer 
Planet Entry Probes. AIAA Paper No. 78-189, J 2 ui. 1978. 

27. Tiwari, S. N.; and Subramauiian, S. V.> Significauice of Radiation 
Models in Investigating the Flow Phenomena Around a Jovian Entry 
Body. Technical Report NASl-14193-26, Jam. 1978, School of 
Engineering, Old Dominion University, Norfolk, Virginia. Also, 

AIAA Paper N6. 78-188, Jan. 1978. 

s. 

28. Vincenti, W. G.; auid Kruger, C. H.i Introduction to Physical Gas 
Dynamics, John Wiley and Sons, 1965. 

29. Spairrow, E. H.; and Cess, R. D.; Radiation Heat Transfer , Brooks/ 

Cole Publishing Co., Belmont, California, 1966. 

30. Penner, S. S.; amd Olfe, D. B. : Radiation amd Reentry , Academic 
Press, 1968. 

31. Hayes, W. D.;and Probstein, R. F.J Hypersonic Flow Theory , 

Academic Press, New York, 1959. 

32. Cheng, H. K. : The Blunt-Body Problem in Hypersonic Flow at Low 

Reynolds Nimber. Inst. Aerospace Science Paper, No. 63-92, 1963. 

33. ShidlovsJcy, V. P.: Introduction to Ratrefied Gases , New York, 

Americam Elsevier Publishing Compauiy, Inc., 1967. 

34. Probstein, R. F.;and Pan, Y. S. : Shock Structure and Leading 

Edge Problem. Rarefied Gas Dynaunics , Vol. 2 CLaurmamn, J. A. Ed.), 
New York, Academic Press, 1963. 

35. Pan, Y. S.;and Probstein, R. F. : Rarefied-Flow Transition at a Leading 

Edge. Proceedings of International Symposium on Fundamental Phenomena 
in Hypersonic Flow, Cornell Univ. Press, 1966, pp. 259-306. 

36. Rott, N.; and Lenard, M. : The Effect of Slip, Particularly for 

Highly Cooled Walls. Journal Aerospace Science, Vol. 29, 1969, 
pp. 591-595. 

37. Anon: The Planet Jupiter, NASA SP-8069, Dec. 1971. 

38. Hoss, J. N. ; Anderson, E. C.; and Bolz, C. W. : Aerothermal 

Environment for Jovian Entry Probes. AIAA Paper No. 76-469, 

June 1976. 

39. Sutton, K. ; Jones, J. J.f amd Powell, R. W. ! Effects of Atmo- 
spheric Structure on Radiative Heating for Jupiter Entry Probe. 

AIAA Paper, No. 78-185, Jan. 1978. 


169 



40 


, Shapiro, A. H. : The Dynamics ^uld Thermodynaipics of Compresalble 
Fluid Flow , Vol. I, The Ronald Press, New York, 1953. 

41. Liepoumn, H. W.; and Roshko, A,; Element of Gaedynami cs , John 
Wiley and Sons, 1957. 

42. Tiwari, S. N,; and Szema, K. Y. : Radiation Induced Precursor 

Plow Field Ahead of a Jovian Entry Body. NASA CR-145221, 1977. 

43. Tiweuri, S. N.; 2 md Szema, K. Y.; Radiation-^ Induced Precursor 
Flow Field Ahead of a Jovieui Entry Body, AIAA Paper No. 77-768, 
June 1977. 

44. Blottner, F. G. i Finite Difference Methods of Solution of the 
Boundary-Layer Equation, AIAA Jouimal , Vol. 8, No. 2, Feb. 1970, 
pp. 193-205. 

45. Tiwari, S. N.;and Szema, K. Y.l Influence of Precursor Heating 
on Viscous Flow Around a Jovian Entry Body. NASA CR-3174, 1979. 

46. Tiwari, S. N.;and Szema, K. Y.; Influence of Precursor Heating 

on Viscous Flow Around a Jovian Entry Body. AIAA Paper No. 78-190, 
Jan. 1978. 

47. Zoby, E. V.; Gnoffo, P. A.; cuid Graves; R. A.: Correlations for 

Determining Thermodyneunic Properites of Hydrogen-Helitm Gas Mix- 
tures at Temperatures from 7,000 to 35,000 K. NASA TN D-8262, 

Aug. 1976. 

48. Street, R. E. ; Problem of Slip Flow in Aerodynamics, NASA RH 
S-7A30, 1957. 

49. Street, R. E.; A Study of Boundairy Conduction in Slip-Flow 
Aerodynamics, Rarefied Gas Dynamics (Devienne, F. M. Ed.), 

London Pergamon Press, 1960, pp. 276-292. 

50. Davis, R. T.; euid Fluge-Lotz, I.; Second Order Boundary-Layer 
Effect in I^ersonic Flow Past Axisymmetric Blunt Bodies. 

Journal of Fluid Mechauiics , Vol. 20, May 1964, pp. 593-623. 

51. Scott, C. D.; Reacting Shock Layer with Slip and Catalytic 
Boundary Condition, AIAA Jotirnal , Vol. 13, Oct, 1975, 

pp. 1271-1278. 

52. Chow, R. R.; and Ting, L. J.l Aerospace Science . Vol, 28, 1961, 
pp. 428-430. 

53. VauiDyke, M. D.t Perturbation Method in Fluid Mechanics , New York 
Academic Press, 1964. 


170 



54 


. Probstein, R. F.jauid Kemp, H. N. , Jr.; Viscous Aerodynamic 
Characteristic in Hypersonic Rarefied Gas Flow, Aero/Space 
Science, Vol. 27, 1960, pp. 174-193. 

55. Ditchbum, R. W. ; Absorption of Ultra-Violet Radiation by the 
Atmospheric Gases . Proceedings of the Royal Society, Series A, 
Mathematical and Physical Sciences , No. 1205, Vol. 236, Aug. 1956, 

pp. 216-226. 

56. Babu, S. G. : Approximate Thermochemical Tables of Some C-H 

and C-H-0 Species. NASA CR-2187, March 1973. 

57. McBridge, B. J.} Heimel, S.; Ehlers, J. G.;and Gordon, S.S 

Thermodynamic Properties to 6000 K for 210 Substances Involving 
the First 18 Element. NASA Report SP-3001, 1963. 

58. Bird, R. B. ; Stewart, W. E. ; and Lightfoot; E. N.: Transport 
Phenomena , John Wiley and Sons, 1960. 

59. Hall, N. A.; Thermodynamics of Fluid Flow , Prentice Hall, Inc., 
1957. 

60. Esch, D. D. ; Pike, R. W. ; Engle, C. D. ; Farmer, R. C. ; and 

Balhoff , J. F. : Stagnation Region Heating of a Phenolic Nylon 

Ablator During Return from Plauietaury Mission. NASA CR -112 026 , 
Sept. 1971. 

61. Leibowitz, L. P. I Measurement of the Structure of eui Ionizing 
Shock Wave in a Hydrogen-Helium Mixture. The Physics of Fluids , 
Vol. 16, No. 1, Jan. 1973, pp. 59-68. 

62. Tiweuri, S. N.; and Szema, K. Y. ! Effects of Precursor Heating on 
Chemical and Radiative Nonequilibrium Viscous Flow Around a Jovian 
Entry Body. AIAA Paper No. 78-907, May 1978. 

63. Tiwari, S. N.;and Szema, K. Y.I Effects of Precursor Heating on 
Chemical and Radiative Nonequilibrium Viscous Flow Around a Jovian 
Entry Body. NASA CR-158132, 1978. 

64. Blottner, F. G.: Viscous Shock Layer at the Stagnation Point with 

Nonequilibrium Air Chemistry. AIAA Journal , Vol. 7, No. 12, 

Dec. 1969, pp. 2281-2288. 

65. Appleton, J. P.; and Bray, K. N. C.: The Conservation Equations 

for a Nonequilibriiam Plasma. Journal of Fluid Mechcuiics , Vol. 20, 
Part 4, Dec. 1964, pp. 659-672. 

66. Kennet, H.; and St rack, S. L.; Stagnation Point Radiative Heat 
Transfer, ARS J. , Vol. 31, No. 3, 1961, pp. 370-372. 

67. Hoshizaki, H.; and Lasher, L. E.t Convective eind Radiative Heat 
Transfer to an Ablating Body. AIAA Journal , Vol. 6, No. 8, 1968, 
pp. 1441-1449, 


171 



68 . 


Oiien, Kuei-Yuan: Application of the S Methods to Spherically 

Symmetric Radiative->Transfer Problems. AIAA Paper Mo. 71»466, 

1971. 

69. Wilson, K. H. ; Evaluation of One-Dimensional Approximations for 
Radiative Transport in Blunt Body Shock Layers. NASA CR-1990, 

1972. 

70. Page, w. A. J Aerodynamic Heating for Probe Vehicles Entering 
the Outer Planets. AAS Paper No. AAS-71-144, 1971. 

71. Lee, P.;and Weissler, G. L. : Absolute Absorption of the 

Continutan, Astrophysical J. , Vol. 115, 1952, pp. 570-571. 

72. Cook, G. R.;and Metzer, P. H. : Photoionization auid Absorption 

Cross Section of and in the Vacuum Ultraviolet Region, 

2 2 • 

J. Opt. Soc. America , Vol. '54, No. 8, Aug. 1964, pp. 968-972. 

73. Samson, J. A. R.; and Cairns, R. B.t Total Absorption Cross 
Section of N^/ and O^ in the Region 550-200 A,. J. Opt. 

Soc . America , Vol. 55, No. 8, Aug. 1965, p. 1035. 

74. McCarty, R. D. : Hydrogen Technology Siirve^i Thermophysical 

Properties. NASA Report SP-3089, 1975- 

75. Gerald, C. F.; Applied Numerical Analysis , Addison Wesley, 1970. 

76. Grose, W. L.;and Nealy, J. E.: Estimates of Noneguilibrium 

Radiation for Venus Entry, AIAA Journal , Vol. 13, No. 4, 

April 1975, pp. 421-424. 


*U.S. GOVERNMENT PRINTING OPFICB: 1979 - 635-004/32 



1. Report No. 2. Government Accession No. 

NASA CR- 3186 

3. Recipient's Catalog No. 

4 . Title and Subtitle 

Effects of Precursor Heating on Radiative 
and Chemically Reacting Viscous Flow 
Around a Jovian Entry Body 

5. Report Date 

October 1979 

6. Performing Organization Code 

7. Author(s) 

S. N. Tiwari and K. Y. Szema 

8. Performing Organization Report No. 

10. Work Unit No. 

9. Performing Organization Name and Address^ 

Old Dominion ; University , Re search Foundation 

P. 0. Box 6369 

Norfolk, Virginia 23508 

11. Contract or Grant No. 
NSG-1492 

13. Type of Report and Period Covered 

Contractor Report 

Feb. 1, 1978-Mar. 31, 1979 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

Langley Technical Monitor; Randolph A, Graves, Jr. 
Final Report 

16. Abstract 


The influence of changes in the precursor region flow properties on the entire 
shock- layer flow phenomena around a Jovian entry body is investigated under physically 
realistic conditions. In the precursor region, the flow is considered to be inviscid 
and the variations in flow properties are determined by employing the small perturba- 
tion technique as well as the thin- layer approximation. The flow in the shock layer 
is assumed to be steady, axisymmetric, and viscous. The analysis is carried out by con- 
sidering both the chemical equilibrium and nonequilibrium composition of the shock- 
layer gas. The effects of transitional range behavior are included in the analysis of 
high altitude entry conditions. 

Realistic thermo-physical and radiation models are used, and results are obtained 
by employing the implicit finite difference technique in the shock layer and an itera- 
tive procedure for the entire shock layer-precursor zone. Results obtained fora 45 
hyperboloid blunt body entering Jupiter's atmosphere at zero angle of attack indicate 
that preheating of the gas significantly increases the static pressure and temperature 
ahead of the shock for entry velocities exceeding 36 km/sec. The nonequilibrium radia- 
tive heating rate to the body is found to be significantly higher than the corresponding 
equilibrium heating. The precursor heating, in general, increases the radiative and 
convective heating to the body, and this increase is slightly higher for the non- 
equilibrium conditions. 


17. Key Words (Suggested by Author(s)) 

Shock- Layer Precursor Heating, 
Nonequi librium Radiative 
Heating, Heating of Entry Bodies 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 34 


19. Security Qassif. (of this report! 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 
180 


22. Price* 


$9.00 


For sale by the National Technical Information Service, Springfield, Virginia 22161 



